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0.0  INTRODUCTION  AND  SUMMARY 

J.  S.  LEE  ASSOCIATES,  INC.  (LAI)  has  shown  in  a  previous  study 
[1],  under  contract  N00014-77-C-0056,  that  an  FM  correlator  is  a  viable 
spectral  line  detector  that  can  be  utilized  for  detection  and  classification 
of  undersea  targets.  The  FM  correlator  was  proposed  as  a  cost-effective 
alternative  tothe  conventional  scheme  which  computes  discrete  Fourier 
transforms  (DFT)  based  on  the  observation  of  a  time  series  over  a  finite 
time  interval.  The  methods  presently  employed  in  spectral  estimation 
for  the  purpose  of  detecting  the  targets'  spectral  band  are  based  on  DFT. 

In  some  applications,  these  methods  are  not  practical  from  the  point 
of  view  of  complexity  and  cost.  For  example,  the  design  of  expendable 
sensor  array  is  usually  constrained  by  a  limited  cost,  and  hence  it  is 
desirable  to  employ  a  scheme  which  is  inherently  simple  and  relatively 
inexpensive. 

The  study  [1]  has  demonstrated  numerically  the  performance  measures 
that  are  quite  acceptable  when  the  FM  correlator  is  designed  under 
suitable  design  parameter  choices.  It  was  observed  that  the  probabilities 
of  detection  can  be  made  to  approach  close  to  unity  for  arbitrary  false 
alarm  rates  when  the  design  parameters  are  chosen  in  a  carefully 
coordinated  manner. 

Having  established  the  concept  validation  of  utilizing  the  FM  correlator 
for  spectral  band  detection,  LAI  has  then  undertaken  the  study  of  designing 
a  total  system  for  effectively  transferring  the  spectral  data  to  a  remote 
processing  facility  via  satellite  relay  channel. 


1 


-  \ 


The  object  of  this  report  is  to  present  design  considerations 
for  the  Spectral  Data  Transfer  System  as  shown  in  an  overall  block  dia¬ 
gram  given  in  Figure  0.1.  Each  subsystem  has  been  analyzed  in  detail, 
with  results  and  recommendation  upon  which  system  designs  can  be  based. 

•  Each  FM  correlator,  employed  as  spectral  line  detector 
makes  a  measurement  in  every  10-15  seconds.  Based  on 
the  resolution  requirement,  these  measurements  are 
quantized  into  approximately  10  bits. 

•  In  view  of  the  low  data  rate  from  the  spectral  detector, 
data  from  a  number  of  the  detectors  will  be  multiplexed 
and  transmitted  in  a  short  burst,  typically  at  1.2  kbps 
or  2.4  kbps. 

«  Scrambled  transmission  is  mandated  in  military  application. 
With  respect  to  the  scrambled  transmission,  a  self¬ 
synchronizing  scheme  is  described.  Self-synchronization 
means  that  synchronization  can  be  recovered  without 
depending  on  a  separate  channel  transmitting  the  sync 
information.  This  scheme,  believed  to  be  new,  is  ex¬ 
plained  in  detail  for  those  without  background  in  the 
theory  of  finite  fields. 

•  Error  correcting  codes  are  briefly  reviewed.  It  is 
concluded  that  either  short  block  codes  with  simple 
threshold  decoding  or  convolutional  codes  with  short 
constraint  length  using  Viterbi  decoding  is  adequate 
for  the  Spectral  Oata  Transfer  System. 
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Figure  0.1  SPECTRAL  DATA  TRANSFER  SYSTEM 


•  Platform  instability  which  severely  restricts  the 
antenna  gain  within  0-3  dB,  and  the  prime  power 
limitation  dictate  the  choice  of  satellites  operating 
at  low  carrier  frequency.  In  order  to  minimize  the 
free-space  loss..  Link  calculations  with  performance 
comparison  have  been  analyzed. 

In  this  report  it  is  shown  that  a  Spectral  Data  Transfer  System 
utilizing  FM  correlators  as  line  detectors,  can  be  designed  which  will 
provide  a  reliable  data  transmission  with  dependable  scrambler  synchroni¬ 
zation  scheme.  The  approach  we  have  taken  in  connection  with  data 
scrambler  and  self  synchronization  scheme  merits  complete  understanding 
of  their  methodology,  for  their  engineering  design  approach  can  be  varied 
in  implementation.  For  these  reasons  we  have  provided  numerous  examples 
as  we  discussed  the  design  philosophies. 
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1.0  SPECTRAL  LINE  DETECTOR 

1.1  The  Concept  of  Using  an  FM  Correlator  for  Line  Detection 

Developing  the  capability  to  detect  and  track  signal  energy 
contained  in  selected  narrow  spectral  windows  continues  to  be  a  significant 
challenge  to  the  undersea-warfare  community.  Potential  targets  can  be 
characterized  acoustically  by  spectra  featuring  one  or  more  narrowband 
emissions  of  uncertain  or  random  bandwidth,  whose  center  frequencies 
are  subject  to  doppler  shifting  as  the  targets  move.  Therefore,  detection 
of  such  emissions  and  estimation  of  their  frequencies  as  they  vary  in  time 
(tracking)  allow  both  target  identification  and  localization.  In  many  ways 
the  dynamic  behavior  of  these  emissions  resembles  frequency  modulation. 

In  this  report  the  design  of  systems  to  exploit  this  resemblance  is  presented. 

A  previous  study  [1]  evaluated  the  performance  of  an  FM  correlator 
in  detecting  the  presence  of  spectral  lines.  The  name  "FM  correlator"  refers 
to  the  signal  processing  configuration  diagrammed  in  Figure  1.1,  in  which 
the  outputs  of  two  sensors  are  first  passed  through  bandpass  filters  (BPF), 
then  demodulated  as  if  they  were  frequency-modulated  (FM)  signals,  then 
correlated.  The  correlation  operation  is  performed  by  multiplying  the  FM 
detector  (FMD)  outputs  and  then  integrating  the  product  by  means  of  a 
lowpass  filter  (LPF).  In  the  analysis  of  this  configuration  the  bandpass 
filter  outputs  were  modelled  by  the  frequency-modulated  waveforms 

si(t)  +  n-(t)  =  Ai  cos[2irfQt  +  (t)]  +  n.. (t) ,  i  =  1,  2 

1  (1-1) 

in  which  the  amplitudes  A.  were  assumed  to  be  constant  and  the  noise  terms 
n.(t)  assumed  to  be  from  independent  stationary,  zero-mean,  narrowband 
Gaussian  random  processes.  The  angle  functions  (t)  were  assumed  to  be 
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Figure  1.1  FM  CORRELATOR  CONFIGURATION 


given  by 


V(t)=|  d^iU)  +  DU)],  i  =  l,2 

(l-2a) 

where 

*  J 

n^t)  <=  m^(t-At)  =  m(t-At) 

( 1 -2b) 

and  the  modulation  m(t)  is  from  a  zero-mean,  stationary,  Gaussian  random 
process.  D(t)  represents  a  slowly  varying  doppler  frequency. shift. 

The  previous  study  indicated  that  successful  detection  based  on 
the  output  of  ..the  FM  correlator  can  be  accomplished  for  certain  values  of 
the  system  parameters.  In  the  following  pages,  further  details  are  given 
concerning  the  system  models  and  the  analytical  results  which  apply  to 
them.  Then,  starting  with  Section.  1.2,  we  begin  to  specify  how  the  FM 
correlator  may  be  used  to  measure  line  frequency  remotely. for  transmission 
to  a  central  data  processing  site. 

1.1.1  FM  Detector  Model 

FM  detectors  (FMD),  or  frequency-to-amplitude  conversion  circuits, 
are  generally  nonlinear  but  so  designed  as  to  Drovide  a  characteristic  that 
is  nearly  linear  for  an  interval  of  frequencies  around  some  center  frequency 
fQ,  as  illustrated  in  Figure  1.2.  Many  FMD  circuits  exist;  an  example  [2] 
of  an  FMD  suitable  for  integrated  circuit  implementation  is  shown  in  Figure 
1.3.  This  circuit  features  less  than  1.5%  departure  from  linearity  for 
frequencies  within  one  perxent  of  the  center  frequency. 

Over  the  linear  region  of  the  characteristic  we  may  express  the 
FMD  outputs  by  the  general  expression 

zi ( t)  -  zQ  =  k. [u>i (t)  -  Uq].  ( 1  -3a ) 

Here  we  assume  the  voltage  offset  Zq  =  0  and  take  the  instantaneous 
angular  frequencies  to  be  derivatives  of  the  phases  shown  in  equations 
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Figure  1.3  SCHEMATIC  OF  INTEGRATED  CIRCUIT 
FM  DETECTOR 


1 


(1-1)  and  (1-2);  that  is. 

“i(t)  ■  "0  * 

(T -3b) 

Then 

zi(t)  =  K- CD  +  m. (t)]. 

(1-4) 

In  this  expression  we  take  D  =  constant  (over  some  interval  of  time)  since 
it  is  slowly  varying. 

1.1.2  Signal  and  Noise  Spectra 

The  noise  components  of  the  received  waveforms  are  assumed  to  have 
the  autocorrelation  functions 

Rn  (T)  =  N.  exp|-(irWNT)2j  COSiOqT  (1-5) 

which  correspond  to  the  Gaussian-shaped  frequency  spectrum 

Sn  (f)  =  — exp{-(f-fQ)2/W2}  ,  f  >  0.  (1-6) 

i  W^/ir 

These  spectra  are  chosen  to  represent  the  class  of  narrowband  spectra, 
and  it  is  assumed  that  WN«fQ.  The  parameter  WN  is  the  one-sided  4.34  dB 
bandwidth;  that  is, 

S  (f±Wj  =  N./e  =  N.  -  4.34  dB.  (1-7) 

n.j  n  i  1 

The  same  spectral  shape  is  chosen  to  represent  the  lowpass  random 
modulation  m(t).  Its  autocorrelation  function  is  assumed  to  be 

'  Pm  «p{-('V>2}-  l’-8) 

1.1.3  The  Correlator  Output 

The  previous  analysis  showed  that  at  the  end  of  a  T-second  observation 
interval,  the  output  of  the  FM  correlator  shown  in  Figure  1.1  can  be 
wri tten 

y(T)  =  u(T)  +  n(T)  (1-9) 

where  y(T)  is  a  mean  value,  given  by 
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-PNR  -PNR  yJ 

R|n(At)  (1-e  -1)(1-e  2)J^  dt  h(t) K2  (1-10) 


where  h(t)  is  the  impulse  response  function  of  the  lowpass  filter  and 
the  CNR^  are  the  carrier-to-noise  power  ratios  at  the  FMD  inputs.  For 
high  CNR.,  then. 


U(T) 


=  const  *  P_ 


m 


(1-11) 


D=0 
I  At=0 

that  is,  the  mean  value  of  the  correlator  output  is  proportional  to  the 
frequency  modulation  power.  For  nonzero  doppler  D  constant  over  the 
T-second  interval ,  we  may  extend  the  result  (1-10)  to  obtain 


p(T) 


At=0 


=  const* (Pm  +  D2). 
m 


(1-12) 


1.2  Using  the  Correlator  Output  to  Measure  Frequency 

The  form  of  the  FM  correlator  output, 

y(T)  =  const  (Pm  +  D2)  +  noise,  (1-13) 

suggests  that  it  may  be  used  to  measure  the  doppler  component  D  when 

2 

the  modulation  power  Pm  is  small.  The  output  depends  on  D  ,  however, 
and  at  most  the  absolute  value  |D|  of  the  doppler  can  be  observed. 

A  simple  modification  to  the  correlator  is  useful.  Consider 
Figure  1.4.  If  an  appropriately  scaled  "offset"  is  added  to  the  FMD 
outputs,  then  the  correlator  output  mean  value  becomes 

p(T)  =  const.  [Pm  +  (D+  n)2]  ~  k(D+  fi)2  (1-14) 

for  P„  small.  This  modification  allows  observation  of  D  including  its 
m  3 

sign,  as  illustrated  in  Figure  1.5. 
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Figure  1.4  FM  CORRELATOR  WITH  OFFSETS 
INTRODUCED  AT  FMD  OUTPUTS 


1 .3  Specifying  the  Frequency  Measurement  System 

The  object  of  this  report  is  to  present  design  considerations 
for  a  system  to  utilize  the  FM  correlator,  assumed  to  be  located  remotely 
with  the  sensors,  to  acquire  measurements  or  estimates  of  target  line 
frequency.  The  overall  diagram  of  such  a  system  is  given  in  Figure  1.6. 
Treatment  of  the  system  components  which  follow  the  FM  correlator  begins 
in  the  next  chapter.  In  this  section,  the  basic  parameters  of  the  FM 
correlator  are  specified. 

1-3.1  Required  Passband 

The  bandwidth  of  a  randomly  modulated  FM  signal  has  been  studied 

2 

by  Abramson  [3],  who  defines  a  mean  square  signal  bandwidth  Bs  by 


Jdf  (f-f0)2ss(f) 
fdf  S$(f) 


(1-15) 


where  S$(f)  is  the  spectrum  of  the  signal.  An  alternate  form  of  (1-15) 
is  given  by 


ul>2 ' 


-HK 


s2rs(t) 


T  =  0 


(1-16) 


in  which  R$(t)  is  the  signal  envelope  autocorrelation  function: 


2 

Rs<x)  ■  TexpK  (0)  +  R*  (t)}- 
'  Ym  Tm 


(1-17) 


Tm  Tm 

From  these  expressions  we  find  that  the  bandwidth  of  the  line  is 


\fm 

2tt 


(1-18) 


Because  of  doppler,  however,  there  is  uncertainty  concerning  where 
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Figure  1.6  REMOTE  FREQUENCY-ACQUISITION  SYSTEM 


in  frequency  this  bandwidth  is  located.  Figure  1.7  illustrates  how  the 

signal  spectum  might  appear  as  a  function  of  time.  It  is  evident  that 

(barfing  any  frequency  tracking  procedure)  the  system  input  passband  is 

determined  by  Dmax,  the  maximum  doppler  shift,  when  the  random  frequency 

modulation  power  P  is  small, 
m 

If  it  is  assumed  that  the  relative  velocity  of  the  acoustic  source 
is  limited  to  within  ±30  knots,  then  it  is  true  that 
2irfn 

l°l  <10#  =  Dmax‘  f1"19) 

and  the  passband  required  at  the  inputs  is 

fQ±U,  so  that  WN  =  fQ/100.  (1-20) 


1.3.2  Smoothing  Requirements 

As  discussed  in  later  sections,  the  noise  remaining  in  the 
correlator  output  is  a  factor  in  the  design.  The  noise  power  is  inversely 
proportional  to  the  time-bandwidth  product  W^T,  where  is  the  (noise) 
bandwidth  at  the  BPF  outputs  and  T  is  the  integration  time. 

The  lowpass  filter  is  chosen  to  perform  an  integrating  or  smoothing 
function.  Ideally,  its  impulse  response  would  be 


h(t)  = 


1 , 0  <  t  <  T 
0  elsewhere. 


(1-21) 


A  more  practical  LPF,  however,  would  have  a  response  that  approximates  this 
function.  For  example,  a  single-pole  RC  filter  has  the  impulse  response 

i 


e-t/RC> 


t  >  0 


h(t)  = 


0-22a) 


0  elsewhere. 


16 


V  s 


with  bandwidth 


B-  sr  o-ao 

The  filter  would  be  "dumped"  or  reset  at  the  beginning  of  the  observation 
interval  and  its  output  sampled  T  seconds  later.  It  was  shown  in  the 
previous  study  [1]  that  the  filter  bandwidth- time  product  BT  has  an  optimal 
value  of  0.4. 
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2.0  QUANTIZATION  AND  A/D  PARAMETERS 


2.1  Resolution  Requirement 

From  external  data  it  has  been  determined  that  a  suitable 
step  size  or  quantization  level  for  f(t)  is 

qf  “2000 

or  qD  =  2.qf  =  D^/20,  since  0^  =  2.fQ/100. 


(2-1) 


What  quantization  level  then  is  required  on  the  correlator  output 
to  insure  this  maximum  q^?  If  we  are  constrained  to  use  uniform 
quantization  of  the  output,  how  many  quantization  levels  are  needed-- 
and  thus  how  many  bits  of  A/D?  These  parameters  are  determined  as 
follows. 

Consider  Figure  2.1.  For  the  no- offset  case,  the  maximum  q^  is 
found  from 


Because  the  slope  of  the  characteristic  is  so  small  near  D=0,  the 
correlator  output  needs  to  be  quantized  very  finely  in  order  to  observe 
the  equivalent  input  doppler  shift.  Thus  800  levels  are  needed  to 
represent  20  levels. 

This  effect  is  even  more  pronounced  for  the  offset  case.  Near  D  = 
we  see  that  the  maximum  q  is  found  from 
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Here,  3200  levels  are  needed  to  represent  20  levels! 

The  situation  can  be  summarized  in  the  following  table: 

Table  2.1  Summary  of  Quantization  Parameters _ 


input  range 

#input  levels 

#output  levels 

bit  "cost" 

o<Id/dbJ<  i 
(no  offset) 

20  (5  bits) 

800  (10  bits) 

5  bits 

-1  <D/D  <1 
'  max 

(offset) 

40  (6  bits) 

3200  (12  bits) 

6  bits 

Because  of  the  nonlinearity,  5  or  6  bits  extra  are  required  to  preserve 
the  measurement  resolution  on  the  frequency.  A  nonuniform  quantizing 
scheme  would  be  more  efficient,  but  also  more  complex.  Therefore,  since 
the  data  rate  is  quite  low  (about  one  sample  every  10  or  15  seconds), 
we  are  willing  to  pay  the  bit  "cost"  to  retain  quantization  simplicity. 

The  fact  that  smaller  values  of  doppler  shift  are  much  more 
likely  than  larger  ones  suggests  a  modified  quantization  strategy. 

Suppose  that  we  do  the  following:  restrict  the  uniform  quantization  of 
the  correlator  output  to  a  smaller  interval  which  corresponds  to 

0  £  |t>/aDMX|  <1,  a  5  1  <2-4> 

then,  for  values  outside  that  interval,  simply  allow  the  A/D  to 
indicate  "overflow",  telling  us  that  a  very  high  doppler  has  been 
observed.  What  effect  does  this  approach  have? 

Consider  Figure  2.1  again.  Using  the  definitions  given  in  that  figure. 


we  can  construct  Table  2.2. 


TABLE  2.2  Quantization  Parameters  for  Restricted  Input  Range 


a 

input 

range 

linput 

levels 

% 

output 
ranqe  - 

^output 

levels 

A/D  bits 

bit  cost 

1 

1.0 

20 

1/800 

1.0 

800 

10 

5 

.9 

.9 

18 

II 

.81 

648 

10 

5 

.8 

.8 

16 

II 

.64 

512 

9 

5 

.7 

.7 

14 

II 

.49 

392 

9 

5 

NO  OFFSET 

.6 

.6 

12 

II 

.36 

288 

9 

5 

.5 

.5 

10 

II 

.25 

200 

8 

4 

1 

2.0 

40 

1/3200 

1.0 

3200 

12 

6 

.9 

1.8 

36 

1/400 

.9 

360 

9 

3 

.8 

1.6 

32 

1/200 

.8 

160 

8 

3 

OFFSET 

.7 

1.4 

28 

3/400 

.7 

94 

7 

2 

Dmax 

.6 

1.2 

24 

1/100 

.6 

60 

6 

1 

.5 

1.0 

20 

1/80 

.5 

40 

6 

1 

The  meaning  of  this  table  is  that,  when  the  output  interval  over  which 
A/D  conversion  is  performed,  is  reduced, 

(a)  for  no  offset,  the  quantization  level  does  not  change,  but 
fewer  levels  are  required  because  the  output  range  is  smaller; 

(b)  for  offset,  the  quantization  level  increases  anc[  the  output 
range  is  smaller,  giving  a  two-way  reduction  in  the  number  of  levels 
requi red . 

Thus,  for  example,  by  deciding  to  declare  "overflow"  for  D  greater  than 
80%  of  its  maximum  value,  we  can  reduce  the  number  of  A/D  bits  required 
from  10  or  12  to  9  or  8.  Under  this  approach,  the  offset  scheme  looks 
very  attractive  because  not  only  does  it  require  one  bit  less  (8  bits),  but 
also  allows  us  to  recover  the  sign  of  the  doppler  shift. 


2.2  Correlator  Output  Noise 

Another  factor  in  determining  the  size  of  the  quantization  levels  to 
be  used  at  the  correlator  output  is  the  amount  of  noise  variation  in  the  output 

of  the  correlator.  For  sufficient  integration  time,  however,  the  noise 
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"ripple"  in  the  output  can  be  kept  very  small  compared  to  the  quantization 
level.  The  amount  of  integration  time  required  is  included  in  the 
discussion  in  Appendix  A. 

2.3  Recommended  Quantization  Scheme 

A  recommended  quantization  scheme  is  given  in  detail  in  the 
Appendix  A  (Table  A-2)  and  features  7-bit  quantization  of  the  reduced  input 

range: 

l°l  *  -75  »«,  (2-5) 

The  20  input  cells  then  correspond  to  unequal  groupinqs  of  output  cells. 

Implementation  of  the  correlator  (including  offset)  and  subsequent 
quantization  operations  is  diagrammed  in  Figure  2.2.  In  the  Appendix  it 
is  shown  that  the  noise  in  the  correlator  output  will  be  sufficiently 
small  to  make  a  128-level  quantization  reasonable  if  T  is  such  that 
WnT  products  on  the  order  of  100  are  produced.  This  is  not  a  stringent 
requirement. 

2:4  Detection  Considerations 

The  output  of  the  FM  correlator  can  be  utilized  for  line  detection, 
suitably  defined.  If  the  output  is  tested  as  to  whether  it  exceeds  a 
threshold,  this  operation  is  equivalent  to  testing  the  hyptothesis 
(for  no  offset) 

Hg:  no  frequency  modulation  present.  (2-6) 

Clearly,  this  test  is  not  equivalent  to  what  is  normally  meant  by 
"detection",  in  which  the  hypothesis  to  be  tested  is  "no  signal  present." 
That  is,  there  is  no  equivalence  unless  it  can  be  stated  positively  that 
the  "signal"  by  definition  contains  frequency  modulation.  If  we  are 
willing  to  make  that  statement,  then  the  analysis  proceeds  as  follows. 

Because  of  the  integrating  effect  of  the  lowpass  filter,  we  can 
treat  the  correlator  output  as  a  Gaussian  random  variable.  Therefore,  we 
have  probabilities  of  false  alarm  and  detection  given  by 
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where  n  is  the  detection  theshold,  and  (m^ o~:)  are  the  means  and  variances 
under  the  hypotheses 

Hq:  CNR  =  0  ,  Hy  CNR>0.  (2-8) 

2  2  2  . 

The  variance  is  given  above  as  oy;  m^p;  mg=0;  and  is  given  by 


2  „2„2  n2  1.35  .299 

° 0  ‘  *1*2  max  2*WT  '  ^^2 


(2-9) 


By  selecting  n  to  guarantee  a  specified  false  alarm  probability  we  may 
relate  PQ  and  PpA  by  the  equation 

XD  °1  +  ml  XFA  °0  (2-10) 


or  xr 


=2ttMT. 


(2-11) 


The  notation  xp  corresponds  to  inverse  probability  and  for  the  Gaussian 
distribution  we  have  the  inverse  distriubtion  given  in  Table  2.3. 

TABLE  2.3  Inverse  Gaussian  Distribution 


p 

xp 

.0001 

3.71902 

.001 

3.09023 

.01 

2.32635 

.5 

0 

.9 

-1.28155 

.99 

-2.32635 

This  equation  can  be  solved  for  a  or  WT,  resulting  in  the  values  given 
in  Table  2.4,  in  which  and  Y  h  P^/UttW^)^. 


TABLE  2.4  Required  WT  Products  for  Detection  Based  on 
FMC  Output 


.5  .. 

.9 

.99 

.5 

.9 

.99 

WT=4.40 

4.42 

4.44 

4.90 

4.67 

4.89 

5.85 

5.87 

5.88 

5.85 

6.12 

6.34 

7.04 

7.06 

7.07 

7.04 

7.31 

7.53 

39.99 

40.01 

40.04 

39.99 

40.29 

40.54 

53.12 

53.14 

53.17 

53.12 

53.42 

53.67 

63.93 

63.95 

63.97 

63.93 

64.23 

64.48 

With  W=fQ/100,  we  can  see  that  an  integration  time  T  of  20  sec  for  low 
speed  targets  (6=1/10)  and  of  2  sec  for  high  speed  targets  (6=1/3)  is 
required  for  PpA=.01. 

Power  Detection.  If  detection  of  the  line  is  to  take  place  in 
the  usual  way  by  square- 1  aw/ lowpass  filtering  as  shown  in  Figure  2.3, 
detection  requires  a  certain  value  of  the  product  WTp(CNR)  [4]  as  given 
in  Table  2.5. 


TABLE  2.5  Required  WTn  and  CNR  For  Power  Detection 

CNR  (dB 


SQUARE 


Obviously,  we  can  make  the  required  CNR  for  detection  as  small  as  we 
please  by  manipulating  the  integration  time  Tp.  This  integration  time 
is  independent  of  that  used  in  the  correlator  and  is  inversely  proportional 
to  the  bandwidth  of  the  lowpass  filter. 
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3.0  SCRAMBLER  AND  SYNCHRONIZATION  SCHEME 


In  this  section  we  first  present  the  proposed  data  scrambler 
system  with  a  brief  description  of  the  main  parts  of  the  system.  We 
then  proceed  to  discuss  each  part  in  detail.  In  developing  a  digital 
data  scrambler  and  descrambler  system,  careful  consideration  must  be 
given  to  the  problem  of  synchronization.  A  self-synchronization  scheme, 
believed  to  be  new,  is  thoroughly  analyzed  both  in  its  principle  of 
operation  and  its  implementation. 

3.1  Description  of  the  Proposed  System 

3.1.1  Scrambler 

A  block  diagram  of  the  scrambler  is  shown  in  Figure  3.1.  The 
scrambler  will  perform  the  following  functions: 

(1)  Receive  and  store  data  from  the  line  detectors 

(2)  Scramble  the  data 

(3)  Arrange  the  data  in  proper  frame  format 

(4)  Provide  signals  for  sync  recovery 

K-2  bits  of  data  are  received  from  each  of  the  N  spectral  detectors  and 
stored  in  K-stage  buffers  with  two  leading  stages  containing  zeros.  At 
the  time  of  transmission  the  multiplexer  switches  out  the  stored  data- 
The  contents  of  each  buffer  are  switched  out  continuously  until  the  buffer 
is  empty.  The  multiplexed  data  are  scrambled  by  combining  with  a  Gold 
sequence  which  is  the  bit-by-bit  modulo-2  sum  of  two  pseudorandom 
sequences,  the  X  sequence  and  the  Y  sequence.  The  scrambled  sequence  is 
blocked  into  frames.  When  A=1  and  B=0,  the  output  is  a  frame  sync  bit; 
when  A=0  and  B=l,  the  output  is  an  original  bit  from  the  X  sequence;  and 
at  the  remaining  instance  both  A=0  and  B=0,  the  output  is  the  scrambled 
data  bit.  The  original  Y  sequence  forms  the  other  part  of  the  output. 
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CLOCK  SWITCH 
AT  BIT 
RATE 


Figure  3.1  DATA  SCRAMBLER 


T 


Both  scrambled  sequence  and  the  Y  sequence  are  sent  to  the  FEC  encoder. 

The  system  is  initialized  by  a  preamble  which  consists  of  a  burst 
of  the  carrier  for  carrier  sync,  a  zero-one  or  all -one. pattern  for  bit 
sync,  and  a  Barker  sequence  or  a  Neuman-Hofman  sequence  for  initialization 
of  frame  sync.  The  circuit  will  be  described  in  Section  3.4. 

3.1.2  Receiver  and  Data  Descrambler 

A  block  diagram  of  the  receiver  and  the  data  descrambler  is  shown 
in  Figure  3.2.  At  the  end  of  the  preamble,  the  receiver  establishes  the 
initial  clock  and  the  initial  sync  for  the  de-interleaver  and  the  pseudo¬ 
random  sequence  generators.  The  de-interleaver  separates  the  Y  sequence 
and  the  scrambled  data  sequence  and  the  frame  sync  circuits  separate 
the  frame  sync  bits  and  the  sampled  bits  of  the  X  sequence  from  the 
scrambled  data  sequence.  Two  pseudorandom  sequence  generators  identical 
with  the  corresponding  ones  in  the  transmitter  and  with  pre-determined 
initial  states  are  activated  to  form  the  Gold  sequence  which  combines 
with  the  scrambled  data  to  recover  the  original  spectral  line  data.  The 
.  received  Y  sequence  and  the  samples  of  the  X  sequence  are  used  for  the 
purpose  of  recovering  the  synchronization  when  sync  loss  is  detected. 

3.1.3  Gold  Sequence  and  Sync  Recovery 

Figure  3.3  shows  two  local  PN  Generators  which, generate  two  pseudo¬ 
random  sequences.  Gold  sequences  are  generated  by  bit-by-bit  modulo-2 
addition  of  the  two  pseudorandom  sequences.  Sync  recovery  for  PN  Generator 
#2  is  accomplished  by  direct  bit-by-bit  comparison  between  the  received 
Y  sequence  and  the  contents  of  the  generator.  Sync  recovery  for  PN  Generator 
#1  is  accomplished  by  a  self-synchronization  circuit  with  the  aid  of  the 
received  samples  of  the  X  sequence.  The  self-synchronization  circuit  will 
be  described  and  analyzed  in  detail  in  Section  3.3. 
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Figure  3.3  GOLD  SEQUENCE  GENERATION  AND  SYNC  RECOVERY 


3.2  Data  Scrambler 

3.2.1  Pseudorandom  Sequences  (A  Tutorial  Discussion) 

In  principle,  data  scrambling  is  achieved  by  changing  the  data 
sequence  "randomly"  before  transmission.  At  the  receiver  the  scrambled 
sequence  is  "changed  back  to  the  original  data  sequence."  The  problem 
lies  in  the  two  phrases:  "randomness"  and  "changing  back".  If  it  is 
completely  random,  the  receiver  will  have  no  way  to  change  back.  On  the 
other  hand  if  the  receiver  knows  how  to  change  back  it  cannot  be  completely 
random.  Consider  the  following  sequences. 

Data  Sequence  .  .  .11001010010101  .  .  .  . 

"Random"  Sequence  .  .  .10100001011010.  ..  . 

Transmitted  Sequence  .  .  .01101011001111  .  .  .  . 

The  transmitted  sequence,  or  the  scrambled  sequence,  is  the  bit-by-bit 
modulo-2  sum  of  the  data  sequence  and  the  "random"  sequence.  At  the 
receiver,  an  identical  "random"  sequence  is  added  to  the  transmitted 
sequence,  yielding  the  original  data  sequence. 

Transmitted  Sequence  .  .  .01101011001111  .  .  .  . 

"Random"  Sequence  .  .  .Iul00001011010.  .  .  . 

Data  Sequence  .  .  .11001010010101  .  .  .  . 

This  simple  illustration  reveals  two  fundamental  requirements  on  the 
"random"  sequence:  it  must  be  reproducible  at  the  receiver  and  it  must 
be  reproduced  in  synchronism  with  the  sequence  which  scrambled  the 
data  sequence.  These  two  fundamental  requirements  make  it  virtually 
impossible  to  use  a  completely  "random"  sequence.  What  is  required,  in 
practice,  is  a  sequence  that  has  sufficient  "randomness"  to  be  unrecoonizable 
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to  the  unintended  observers  and  yet  is  deterministic  making  it  relatively 
easy  to  generate  and  synchronize  in  the  receiver. 

The  most  important  method  of  generating  binary  sequences  is  by 
means  of  a  shift  register  with  feedback  connections.  An  n-stage  shift 
register  is  a  device  consisting  of  n  consecutive  binary  storage  elements. 

The  contents,  either  "0"  or  "1",  of  each  element  can  be  shifted  to  the 
next  position  down  on  receipt  of  a  regular  clock  pulse.  The  feedback 
will  be  a  "1"  or  a  "0"  and  may  be  computed  as  a  logical  function  of  the 
contents  of  the  shift  register  stages.  In  a  linear  feedback  shift 
register,  which  is  the  most  useful  and  practical  realization  of  a  shift- 
register  sequence  generator,  the  feedback  function  is  a  modulo-2  sum  of 
two  or  more  of  the  shift  register  stages.  An  example  of  a  5-stage  linear 
feedback  shift -register  sequence  generator  is  shown  in  Figure  3.4. 

Five  arbitrary  binary  digits,  not  all  zeros,  are  loaded  into  the  shift 
register  as  the  initial  conditions.  At  each  clock  time,  the  content  of  each 
shift  register  stage  shifts  to  the  next  stage;  the  content  of  the 
last  stage  and  the  content  of  the  second  stage  are  modulo-2  added;  and  the  : 
is  then  fed  into  the  first  stage.  At  the  same  time  the  content  of  the 
last  stage  is  taken  out  as  the  output  sequence. 

For  a  shift-register  sequence  generator  of  n  stages,  the  output 
sequence  will  be  always  periodic  because,  whatever  the  initial  condition, 
after  a  number  of  clock  pulses  the  inital  condition  must  eventually  be 
reproduced.  Since  the  maximum  number  of  different  combinations  of  n 
binary  digits  is  2n,  the  period  of  the  output  sequence  cannot  exceed  2n. 

In  the  linear  feedback  shift  register,  the  state  of  all  zeros  in  the  shift 
register  will  produce  a  "0"  as  feedback;  hence  the  state  of  the  shift 
register  will  be  all  zeros  forever.  The  "all-zeros"  state  is  a  degenerate 


state;  it  cannot  appear  in  the  shift  register  if  the  initial  state  is 
not  "all  zeros";  therefore  the  maximum  number  of  possible  states  is  2n-l. 
An  output  sequence  whose  period  is  2n-l  is  called  a  "maximum  length 
sequence"  or  "m-sequence"  for  short.  The  quantity  L=2n-1  is  called  the 
length  of  the  sequence. 

Maximum-length  sequences  are  also  referred  as  "pseudorandom 
sequences"  or  "pseudonoise  sequences"  or  PN  sequences  for  short.  They 
possess  the  following  "randomness"  properties. 

(1)  The  Balance  Property:  In  a  complete  period  of  the  sequence 
the  number  of  "1"  differs  from  the  number  of  "0"  by  at  most  1. 

(2)  The  Run  Property:  There  are  runs  of  consecutive 

"1"  and  "0”,  a  half  of  the  runs  are  of  length  1 ,  1  of  the  runs  are  of 
length  2,  ^  of  the  runs  are  of  length  3,  etc. 

(3)  The  Correlation  Property:  If  a  complete  sequence  is  compared, 
bit  by  bit,  with  any  shift  of  itself  the  number  of  agreements  differs 
from  the  number  of  disagreements  by  at  most  1. 

For  example,  the  sequence  generated  by  the  shift-register  sequence 
generator  in  Figure  3.4  with  1  0  0  0  0  as  the  initial  state  will  be: 

0000101011101100011111001101001. 

5 

The  length  of  the  sequence  is  2-1  =  31.  The  balance  property  is 
satisfied  as  there  are  16  "1"  and  15  "0".  The  run  property  is  satisfied 
as  there  are  — =  16  runs  (0  runs  of  "1"  and  8  runs  of  "0"),  there  are 
8  runs  of  length  1  (4  runs  of  "1"  and  4  runs  of  "0"),  there  are  4  runs 
of  length  2  (2  runs  of  "1"  and  2  runs  of  "0"),  there  are  2  runs  of  length 
3  (1  run  of  "1"  and  1  run  of  "0");  and  finally  there  is  1  run  of  "0"  of 
length  4  and  1  run  of  "1"  of  length  5.  Finally  the  correlation  property 
is  satisfied  as  follows: 
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0  0  0  0  1  0  1  0  1  1 
1  0  0  0  0  1  0  1  0  1 

a  a  a  a 

d  ddddd 


1  0  1  1  0  0  0  1  1  1 

1  1  0  1  1  0  0  0  1  1 

a  a  a  a  a  a 

d  d  d  d 


1  1  0  0  1  1  0  1  0  0  1 
1110  0  110  10  0 

a  a  a  a  a 

d  d  d  d  d  d 


complete  sequence 
shifted  sequence 

15  agreements 

16  disagreements. 


The  comparison  between  a  complete  sequence  with  a  shifted  sequence  shows 
15  agreements  and  16  disagreements. 

The  feedback  connections  of  a  LFSR  (abbreviation  of  Linear  Feedback 
Shift  Register)  of n  stages  can  be  made  into  a  cne-to-one  correspondence 
with  a  polynomial  of  degree  n  with  either  "0"  or  "1"  as  coefficients, 
where  a  coefficient  of  1  represents  a  connection  (including  two  ends) 
and  0  represents  no  connection.  Consider  the  5-stage  sequence  generator 
in  Figure  3.4  and  redraw  it  in  Figure  3.5.  The  corresponding  polynomial 
is 

P(x)  =  lx5  +  Ox4  +  Ox3  +  lx2  +  Ox1  +  lx° 


arranged  in  an  ascending  order  from  left  to  right.  The  polynomial  which 
corresponds  to  the  feedback  connections  of  a  maximum  length  sequence 
generator  or  PN  generator  is  called  the  "primitive  polynomial."  Each 
primitive  polynomial  of  degree  n  will  generate  a  unique  pseudorandom 
sequence  of  length  L  =  2n-l. 

The  next  questions  are:  how  many  primitive  polynomials  are  there 
of  degree  n  and  how  can  we  find  the  primitive  polynomials? 

The  answer  to  the  first  question  is:  there  are  exactly  #(n)/n  primitive 
polynomials  of  degree  n,  where  $(n)  is  the  number  of  integers  larger  than 
0,  less  than  L  =  2n-l,  and  relatively  prime  to  L  (including  1). 
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For  example,  for  n=4,  L=2n-1  =  15,  the  integers 
1,  2,  4,  7,  8,  11,  13,  14 
are  relatively  prime  to  15;  hence 
4>(n)  =  8 

and  the  number  of  different  primitive  polynomials  is 

<fr(n)  _  8  =  2 
n  4 

5 

For  n=5,  since  L=2  -1=31  is  a  prime  number  itself,  all  integers  less 
than  31  are  prime  to  31;  therefore 
$(n)  =30 

and  the  number  of  different  primitive  polynomials  is 

4>(n)  _  30  _  g 
n  “  5  "  6' 

A  list  of  the  number  of  possible  pseudorandom  sequences  for  shift 
registers  having  up  to  21  stages  is  given  in  Table  3.1. 

To  find  the  primitive  polynomials  of  degree  n  turns  out  to  be  a 
very  difficult  problem.  Fortunately,  tables  ot  irreducible  polynomials 
has  been  published  up  to  degree  34,  among  which  the  primitive  polynomials 
are  labeled.  A  well  known  source  is  in  Peterson’s  book  [5,  Appendix  C, 
pp.  472-492].  We  shall  point  out  the  notations  used  in  that  reference  in 
order  to  be  able  to  make  use  of  the  table.  Take  an  entry: 

Degree  5  1  45E  3  75G  5  67H. 

The  numbers  45,  75,  and  67  represent  the  polynomials  in  octal 
representation: 

0  000  2  010  4  100  6  110 
1  001  3  011  5  101  7  111. 

Thus,  45  represents  100101 
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NUMBER  OF 
SHIFT  REGISTER 
STAGES  n 


SEOUENCE 
LENGTH 
L  =  2n-l 


TOTAL  NUMBER 
OF  SEQUENCES 


2 

3 

1 

3 

7 

2 

4 

15 

2 

5 

31 

6 

6 

63 

6 

7 

127 

18 

8 

255 

16 

9 

511 

48 

10 

1023 

60 

11 

2047 

176 

12 

4095 

144 

13 

8191 

630 

14 

16383 

756 

15 

32767 

1800 

16 

65535 

2048 

17 

131071 

7710 

18 

262143 

7716 

19 

534287 

27594 

20 

1048575 

24000 

21 


2097151 


84672 


which  in  turns  represents  the  polynomial 

lx5  +  Ox4  +  Ox3  +  lx2  +  Ox  +  1 
or 


The  number  75  represents  111101 

5  4  3  2 

or  the  polynomial  x  +  x  +  x  +  x  +1 

and  67  represents  110111 

5  4  2 

or  the  polynomial  x  +x  +x  +x+l. 

The  letter  after  the  numbers  gives  the  following  information 
about  the  polynomial:  E,  F,  G,  H  indicate  that  the  polynomial  is 
primitive,  while  A,  B,  C,  D  indicate  not  primitive. 

The  number  before  45,  75,  67,  that  is  1,3,  5,  has  the  following 

meaning.  If  a=a*  is  a  root  of  the  polynomial  represented  by  45,  then 

3  5 

a  is  a  root  of  the  polynomial  75  and  a  is  a  root  of  polynomial  67;  or 

5  2 

a  is  a  root  of  x  +  x  +  1 . 

3  5  4  3  2 

a  is  a  root  ofx  +  x  +  x  +  x  +1, 

5  .  **5J4.2J. 

and  a  is  a  root  ofx  +x  +  x  +x  +  l. 

For  degree  5,  three  primitive  polynomials  are  listed: 

45  -*-►  1  00101  x5  +  x2  +  1 

75*-*  1  1  1  1  0  1  «*-*■  x5  +  x4  +  x3  +  x2  +  1 

67 -*-*-1  1  01  1  1  «-►  x5  +  x4  +  x2  +  x  +  1 . 

They  will  generate  three  different  pseudorandom  sequences  of  length 

L=31.  Each  polynomial  corresponds  to  a  PN  generator.  We  observe  that 

in  all  three  generators  the  content  of  the  last  stage  always  feeds 

back  to  the  first  stage  reflecting  the  fact  that  the  coefficients  of  the 

highest  degree  term  (the  x^  term)  and  the  lowest  degree  term  (the  x®  term 

or  the  constant  term)  always  equal  1.  In  the  first  generator  the  content  of  the 

second  stage  also  feeds  back;  in  the  second  generator  the  contents  of  the 
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second,  third,  and  fourth  stages  feedback;  and  the  contents  of  first, 

second,  and  fourth  stages  feedback  in  the  third  generator.  The  first 

generator  requires  one  tap  while  the  other  two  generators  require 

5  2 

three  taps.  The  generator  corresponding  to  the  polynomial  x  +  x  +  1 
is  called  "the  minimum  tap  pseudorandom-sequence  generator."  In  the 

Table  of  the  Irreducible  Polynomials  the  primitive  polynomial  corresponding 
to  the  minimum-tap  generator  is  always  listed  first. 

According  to  Table  3,1  there  arc-  six  distinct  pseudorandom  sequences, 
yet  there  are  only  three  primitive  polynomials  in  the  Table  of 
Irreducible  Polynomials.  The  three  missing  polynomials  can  be  obtained 
as  follows.  Let 

P1(x)  =  x5  +  x2  +  1. 

Replace  x  by  x~* 

M*’1) =  +  “V  +  l- 

1  x3  x 

-1  5 

Multiplying  P^x  )  by  x  ,  we  obtain 

x^tx-1)  =  1  +  X3  +  x5  =  x5  +  X3  +  1. 

5  3 

The  new  polynomial  x  +  x  +  1  is  called  "the  reciprocal  polynomial" 
of  Pj(x),  and  is  represented  by  the  notation  P^(x).  Thus,  the  reciprocal 
polynomial  of  a  polynomial  Pj(x)  is  defined  as 

p;u)  =  x^u'1). 

It  can  be  shown  that  the  reciprocal  polynomial  of  a  primitive  polynomial 
is  also  primitive. 

To  find  the  reciprocal  polynomials  is,  in  practice,  very  easy 
by  making  the  following  observation: 


10  0  10) 


10  10  0) 


!  I 


v«>  •  *5  ♦  ♦  > 

P?{»)  *  »5  <  »4  ♦  x3  ♦  X2  ♦  1  —  111101 
P3(x)  •  x5  ♦  »4  ♦  *2  ♦  *  ♦  1  —  110111 


"  p;i.)  •  x5  .  x3  .  i 

101111  ~  P*(x)  •  x5  «  x3  «  1*  «  x  *  1 

1110  11  —  9\M  -  x5  *  x*  *  l3  «  x  ♦  1  . 


The  procedure  is:  (1)  write  down  the  binary  sequence  representation  of 

the  given  polynomial,  (2)  read  the  binary  sequence  backwards,  and  (3)  chanqe 
the  sequence  to  the  reciprocal  polynomial. 

We  now  obtain  all  6  pseudorandom  sequence  generators  of 

sequence  length  L=31: 

Pj(x)  =  x5  +  x2  +  1 

P2(x)  =  X5  +  X4  +  X3  +  X2  +  1 

P3(x)  =  X5  +  X4  +  X2  +  X  +  1 

P4(x)  =  x5  +  X3  +  1 

P5(x)  =  x5  +  X3  +  X2  +  X  +  1 

P6(x)  =  X5  +  X4  +  X3  +  X  +  1. 


The  polynomials  P-j (x)  and  P4(x)  are  reciprocal  to  each  other;  so  are  P2(x) 
and  P5(x) ,  and  P3(x)  .and  Pg(x). 

Take  an  example.  For  n=ll  ,  we  obtain  from  the  Table  of 
Irreducible  Polynomials 
(1)  1  4005E 

4005  ■*-*  1  000000001  01 
Pj(x)  =  x11  +  X2  +  1  . 

Read  the  sequence  backwards 


101  000000001-^  5001 
and  P?  (x)  =  x11  +  x^  +  1. 


Read  the  sequence  backwards 

10100100100  W  5111 
and  Pg  (x)  =  x11  +  x9  +  x6  +  x3  +  1. 

(3)  13  4143F 

4143  ^10000110001! 

P13(x)  =  x11  +  X6  +  X5  +  X  +  1 
Read  the  sequence  backwards 

1  1  0  0  0  1  1  0  0  0  0  1  «  6141 
and  P*3(x)  =  x11  +  x10  +  x6  +  x5  +  1. 

The  pseudorandom  sequence  generator  can  be  implemented  with 
flip-flops  and  logic  gates.  The  implementation  of  the  generator  in 
Figure  3.4  is  shown  in  Figure  3.6.  The  flip-flops  used  in  this 
implementation  are  JK  flip-flops  with  additional  inputs  of  PRESET  or 
direct  SET  ( SQ)  and  CLEAR  or  direct  RESET  (Rp).  These  inputs,  which 
function  exactly  as  the  S  and  R  inputs  of  the  SR  flip-flops  do,  are 
used  to  establish  the  initial  state  of  the  shift  register.  After  RESET 
clears  the  register,  the  LOAD  line  is  activated  and  the  initial  states 
D-j ,  D2,  D3,  04,  and  D5  can  be  entered  into  the  register  in  parallel. 

Following  the  loading  operation  the  data  shift  one  position  to  the  right 
with  each  clock  pulse  and  at  the  same  time  the  contents  of  the  second  and 
fifth  stages  feed  back  to  the  first  stage  through  an  EXCLUSIVE  OR.  The 
pseudorandom  sequence  is  read  out  serially  at  the  last  stage. 

We  are  now  ready  to  observe  how  the  maximum-length  sequence  is 
generated.  We  do  this  by  considering  a  series  of  examples. 

Example  1.  P{x)  =  x3  +  x  +  1 

The  shift  register  pseudorandom  sequence  generator  is  shown  in 
Figure  3.7.  Let  F^t,  F2t>  and  Fj^  be  the  contents  of  the  corresponding  stages 
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or, 

t  =  01234  567  =  0 

Fj  10  0  1110  1 

F2  0  1  0  0  1  1  1  0 

F3  10  10  0  111. 

The  pseudorandom  sequence  generated  is 
10  10  0  11. 

Example  2 

P(x)  =  x4  +  x  +  1 

The  PN  generator  is  shown  in  Figure  3.8. 


then  the  values  of  F^ ,  F,,,  F^,  and  F^  at  t  =  1,  2,..., 14  can  be 
calculated  by  matrix  multiplication.  The  results  are  as  follows 
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Figure  3.8  4-STAGE  PSEUDORANDOM 
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t  =  0  1  2  3  4  5  6  7  8  9  10  11  12  13  14 

Fj  1001000111  1  0  1  0  1 

F2  1100100011  1  1  0  1  0 

f3  0110010001  1  1  1  0  1 

F4  1011001000  1  1  1  1  0. 

Example  3 

P(x)  =  x5  +  x2  +  1 

The  5-stage  PN  generator  is  shown  in  Figure  3.9  with  the  shift  register 
contents  and  the  output  pseudorandom  sequence  written  directly  under 
the  generator.  The  mat.  !/,  multiplication  method  which  provides  a  clear 
understanding  of  the  operation  of  a  PN  sequence  generator  as  io  Examples 
1  and  2,  becomes  cumbersome  as  the  number  of  stages  increases. 

We  are  now  ready  to  describe  the  principle  of  a  data  scrambling 
system.  When  a  data  sequence  is  added  onto  a  PN  sequence,  the  resultant 
sequence  becomes  unintelligible  to  an  intercepter.  But  the  original  data 
sequence  can  be  recovered  by  adding  an  identical  PN  Sequence.  A 
conceptual  block  diagram  of  a  data  scrambling  and  descrambling  system  using 
a  pseudorandom  sequence  is  shown  in  Figure  3.10.  The  data  sequence  is 
scrambled  by  a  pseudorandom  sequence  with  a  bit-by-bit  modulo-2  adder. 

At  the  receiving  end  a  synchronized  identical  pseudorandom  sequence  is 
bitrby-bit  modulo-2  added  onto  the  scrambled  sequence  to  recover  the 
original  data  sequence.  The  data  scrambling  system  described  thus  far 
and  a  spread-spectrum  system  are  similar.  The  only  difference  is  that  in 
the  scrambling  system  the  data  sequence  and  pseudorandom  sequence  are  of' 
the  same  rate  while  in  the  spread-spectrum  system  the  pseudorandom  sequence 
is  much  faster  than  the  data  sequence. 
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Example  3 


Figure  3.®  5-STAGE  PSEUDORANDOM  SEQUENCE  GENERATOR 


3.2.2  Gold  Codes 

Pseudorandom  sequences  have  extensive  applications  in  spread- 
spectrum  communication  systems.  In  the  case  of  a  spread-spectrum  system 
the  pseudorandom  sequence  serves  as  an  encoding  mechanism  which,  when 
added  to  the  data  sequence,  results  in  a  wideband  signal  as  shown  in 

Figure  3.11.  If  the  communication  system  is  a  multiple  access  system  each  j 

j 

pseudorandom  sequence  is  used  as  a  code  for  a  particular  user. 

The  usefulness  of  the  pseudorandom  sequences  in  a  spread- spectrum  j 

system  depends  in  large  part  on  their  ideal  autocorrelation  properties.  ] 

! 

One  of  the  randomness  properties  of  the  pseudorandom  sequence  is  the  ] 

i 

correlation  property,  i.e.  if  a  complete  sequence  is  compared,  bit  by 
bit,  with  any  shift  of  itself  the  number  of  agreements  differs  from  the 
number  of  disagreements  by  at  most  1.  From  this  property  it  is  easy  to 
see  that  the  autocorrelation  function  of  a  pseudorandom  sequence  is  of 
the  shape  described  in  Figure  3.12. 

The  cross -correlation  function  between  two  different  pseudorandom 
sequences  of  the  same  length  is,  however,  an  entirely  different  matter.  It 
can  have  high  peaks;  and  to  make  the  matter  worse,  there  is  no  method 
available  to  calculate  the  cross-correlation  function  between  two  pseudo¬ 
random  sequences  except  by  simulation.  For  long  sequences  this  is  not 
possible  even  with  the  fastest  computers.  To  visualize  the  cross-correlation 
problem  and  how  much  is  involved  in  the  simulation,  consider  an  example  as  follows. 

3 

The  sequence  generator  x  +  x  +  1  with  initial  state  1  0  0 
generates 

S]  =  0  0  1  1  1  0  1 

3  2 

and  the  sequence  generator  x  +  x  +  1  with  initial  state  1  1  1  generates 
S2  »  1  1  1  0  0  1  0. 


DATA 

SEQUENCE 


PSEUDORANDOM 

SEQUENCE 


1  0  011  1  01  00  1011  1  1  0  0  1  0  0  0 


TRANSMITTED  0  11  00  101  00  1  01110 

sequence  i  i  n  n  n  I  I  I 


110  1  0  0 


Fioure  3.11  SPREAD-SPECTRUM  WAVEFORM 


The  cross-correlation  is  computed  by  holding  fixed,  shifting  S2 
one  bit  at  a  time,  and  comparing  the  two  sequences  to  obtain  the 
difference  between  the  number  of  agreements  and  the  number  of  disagree¬ 
ments. 


1.  0011101 

+1110010 

110  1111  -v  -5 

2.  0  0  1110  1 

tQ  1  1  1  0  0  1 

0  1  0  0  1  0  0  -*-  +3 

3.  0011101 
+1011100 

1  0  0  0  0  0  1  +  +3 

4.  0011101 
+0101110 

0  110  0  11-*-  -1 

5.  0  0  1110  1 
+  0010111 

0  0  0  1  0  1  0  +  +3 

s.  0  0  1110  1 
+  10  0  10  11 

1010110+  -1 

7.  0011101 

+1100101 

1111000+  -1 
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The  cross-correlation  function  is  plotted  in  Figure  3.13. 

In  Section  3.2.1,  we  have  found  all  six  PN  generators  of  5  stages. 

The  cross-correlations  between  sequences  generated  by  Pj(x)  and  P^{x), 
and  by  P^(x)  and  Pg(x)  are  shown  in  Figure  3.14. 

The  cross-correlation  function  between  two  distinct  pseudorandom 
sequences  is  a  very  important  requirement  in  a  multiple  access  communications 
system  as  a  user's  receiver  might  lock  onto  a  wrong  signal  if  the  cross¬ 
correlation  peak  between  the  user's  pseudorandom  sequence  and  thersequence 
for  the  wrong  signal  is  high  enough  to  exceed  the  threshold  for  synchroni¬ 
zation. 

To  overcome  the  cross-correlation  problem.  Gold  considered  the 
bit-by-bit  modulo-2  sum  of  two  pseudorandom  sequences  of  same  length  but 
generated  by  two  distinct  primitive  polynomials,  Pj(x)  and  p2(x),  as  shown 
in  Figure  3.15.  If  the  length  of  the  two  pseudorandom  sequences  is  2n-l, 
then  the  resultant  sequence  also  repeats  itself  after  2n-l  bits.  Furthermore 
if  one  sequence  is  kept  fixed  and  the  second  sequence  is  shifted  in  time, 
a  different  resultant  sequence  is  generated.  2n-l  different  sequences  can 
be  generated  this  way,  one  for  each  different  time  shift  of  the  second 
sequence.  Joining  the  two  original  pseudorandom  sequences,  altogether 
2n+l  different  sequences  can  be  generated  with  one  pair  of  primitive  polynomials. 
These  sequences  are  sometimes  referred  to  as  Gold  sequences  or  Gold  codes; 
they  are  not  maximal  except  the  two  original  pseudorandom  sequences.  It 
should  be  noted  that  for  n-stage  shift  register  there  are  only  «j>(2n-l)/n 
pseudorandom  sequences  and  yet  a  pair  of  n-stage  shift  registers  can 
generate  2n+l  different  Gold  sequences.  The  increase  in  the  number  of 
available  sequences  is  drastic.  For  instance,  from  Table  3.1,  for  n=10, 

*(2n-l)/n  =  60  and  2n+l  =  1025. 
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As  an  example  of  Gold  sequences  consider  the  case  of  n=3.  There 
are  two  primitive  polynomials 

Pj(x)  =  x3  +  x  +  1  generates  S1  =0011101 
and 

Pj(x)  =  x3  +  x^  +  1  generates  S2  =1110010 
This  pair  generates  nine  Gold  sequences 

S.  =  0011101 

52  *1110010 

53  =1101111 

54  =0100100 

55  =  1  0  0  0  0  0  1 

56  =  0110011 

57  =0001010 
Sg  =  1  0  1  0  1  1  0 
S9  -1111000 

We  have  illustrated  that  the  number  of  available  codes  increases 
from  to  2n+l.  However,  we  have  yet  to  address  to  the  cross¬ 

correlation  problem.  To  solve  this  problem  Gold  proved  the  following  [6]' 
Given  a  pseudorandom  sequence  generated  by  a  polynomial  P.|(x) 

of  degree  n  having  a  as  a  root,  let  P^(x)  be  a  polynomial  with  a*  as 


a  root  where 


if  n  is  odd 


if  n  is  even . 


Then  the  cross-correlation  between  the  Gold  sequences  generated  by  P^(x) 
and  Pt(x)  satisfies  the  following: 


|RT(l,t)U 


+  1  if  n  is  odd 


+  1  if  n  is  even. 


Pj(x)  and  Pt(x)  form  a  preferred  pair. 


>  .  \ 


Figure  3.13  AUTOCORRELATION  FUNCTION  AND  CROSS-CORRELATION 
FUNCTION  BETWEEN  TWO  PSEUDORANDOM  SEQUENCES 
GENERATED  BY  x3  +  x  +  1  and  x3  +  x2  +  1 
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xz  + 

X3  + 
X4  + 


X  +  1 


Figure  3.14  CORRELATION  OF  PSEUDORANDOM  SEQUENCES 
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BIT-BY-BIT 
MODULO- 2 
ADDER 


1 

J 

n-STAGE  PN  GENERATOR 

P2(x) 

Figure  3.15  GOLD  SEQUENCE  GENERATOR 


Suppose  that  n=3, 

3 

Px(x)  =  X  +  X  +  1 
n+1 

?  2 

2^  +  1  =  Zc  +1=5 
P5(x)  =  x3  +  x2  +  1. 

The  absolute  value  of  the  cross-correlation  between  the  Gold  sequences 
3  3  2 

generated  by  x  +  x  +  1  and  x  +  x  +1  does  not  exceed  5. 

We  now  consider  an  example  to  show  how  to  design  a  Gold  sequence 
generator  with  a  preferred  pair  of  primitive  polynomials  of  degree  n. 
Suppose  that  n  is  given: 
n  =  11. 

From  Peterson's  Table  of  Irreducible  Polynomials,  we  can  choose  1  4005E 
as  Pj(x) ,  i .e. 

Pj(x)  =  x11  +  x2  +  1 
n+2 

t  =  2  2  +  1  =  65. 

We  want  to  find  P65(x)  which  has  a65  as  a  root.  But  the  polynomial  with 
65  as  leading  number  is  not  listed  in  the  table.  In  this  case  we  can 
do  the  following. 

(1)  Calculate  2m  x  65,  for  m=l ,2, . . . ,11 .  If  2m  x  65  is  larger 
than  2^-1  =  2047,  divide  2m  x  65  by  2047  and  replace  2m  x  65  by  the 
remainder,  and  then  continue  the  doubling  process. 

G5,  130,  260,  520,  1040,  2080  =  33,  66,  132,  264,  528,  1056,  2112  =  65 

(2)  The  doubling  process  will  end  at  the  ll-th  time;  at  the  12-th 
time  the  number  will  go  back  to  the  original  number,  65. 
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(3)  a65,  a130, . .  .(J056  are  roots  of  the  same  polynomial. 


(4)  One  of  the  eleven  numbers,  33,  is  listed  in  the  table: 

33  7335G. 

(5)  Therefore, 

33' 


P33(x)  -  x11  +  X10  ♦  X9  ♦  X7  ♦  X6  ♦  X4  +  X3  +  x2  +  1. 


,11 


There  are  2  +  1  =  2049  Gold  sequences  of  length  L=2047  with  cross¬ 

correlation 


|R.  .(t) |  <  65. 

* 

Gold  sequences  can  be  used  for  scrambling  just  like  pseudorandom 
sequences.  They  are  especially  attractive  in  systems  which  do  not  require 
long  sequences  but  need  to  accommodate  a  large  number  of  users.,  A 
conceptual  block  diagram  of  a  scrambler  using  a  Gold  sequence  is  shown  in 
Figure  3.16.  The  two  pseudorandom-sequence  generators  are  of  the  same 
number  of  stages  (n  stages,  say).  The  bit-by-bit  modulo-2  sum  of  the 
two  pseudorandom  sequences  is  a  Gold  sequence  which  is  a  (2n-l) -digit 
non-maximal  sequence  (which  means  the  Gold  sequence  is  not  a  maximum- 
length  sequence  although  both  original  pseudorandom  sequences  are). 

3.2.3  Frame  Synchronization 

In  Section  1  of  this  report  we  have  established  the  feasibility 
of  using  FM  correlator  to  measure  the  frequency  of  the  spectral  line. 

The  spectral  line  detector  makes  one  measurement  in  every  10-15  seconds 
and  the  measurement  value  is  quantized  and  encoded  into  a  7-11  bit  word 
for  transmission  to  a  central  location  for  processing.  The  information 
rate  from  a  single  line  detector  is  very  low.  Even  though  a  line  detection 
system  could  combine  a  large  number  of  detectors  the  information  rate 
would  still  be  so  low  that  the  only  viable  transmission  method  would  be  in 
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burst  mode;  i.e.  all  the  spectral  data  are  stored  and  are  being  transmitted 
in  a  short  burst  of  time.  The  typical  transmission  rate  of  a  communication 
link  via  satellite  could  be  either  1200  bps  or  2400  bps.  The  function  of  the 
frame-synchronization  scheme  is  to  separate  the  data  from  different 
detectors  into  blocks  called  frames.  The  start  of  a  frame  is  indicated 
by  one  or  more  bits  periodically  inserted  at  the  beginning  of  each  frame. 

For  the  system  under  discussion,  because  of  the  low  data  rate  a  single  bit 
for  an  alternating  zero-one  pattern  is  sufficient  for  the  receiver  to 
obtain  frame  synchronization.  The  frame  structure  is  shown  in  Figure 
3.17.  It  should  be  noted  that  the  bit  number  is  arbitrarily  chosen  just 
for  illustration  purpose.  Each  frame  contains  the  scrambled  spectral 
data  from  one  FM  line  detector.  At  the  beginning  of  each  frame  one -bit 
from  an  alternating  zero-one  pattern  is  used  for  frame  sync;  the  second 
bit  is  from  the  X-sequence  which  will  be  discussed  in  detail  later;  and 
the  remaining  11  bits  are  the  scrambled  data  from  one  spectral  detector. 

3.2.4  Data  Scrambler  System  and  Block  Diagram 

As  mentioned  in  Section  3.1.1,  the  function  of  the  data  scrambler 
system  is: 

(1)  to  receive  and  store  data  from  the  spectral  detectors 

(2)  to  scramble  the  data 

(3)  to  arrange  the  data  in  proper  frame  format 

(4)  to  provide  signals  for  sync  recovery. 

As  shown  in  the  block  diagram  in  Figure  3.1,  K-2  bits  of  data  are 
received  from  one  of  the  N  spectral  detectorsrand  are  stored  in  a  K-stage 
buffer  with  two  leading  stages  containing  zeros.  At  the  time  of  transmission, 
the  multiplexer  switches  out  the  data  emptying  the  buffers  one  by  one  with 
the  output  clock  running  at  the  same  rate  as  the  bit  rate  clock.  The 
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output  of  the  multiplexer  is  a  continuous  bit  stream  separated  into 
blocks  of  K  bits;  each  block  begins  with  two  zeros  followed  by  K-2  data  bits 
from  one  spectral  detector. 

For  data  scrambling,  a  Gold  sequence  is  chosen  over  a  single 
pseudorandom  sequence.  A  Gold  sequence  possesses  two  advantages  over 
a  single  pseudorandom  sequence:  namely,  the  bounded  cross-correlations 
and  the  large  number  of  available  sequences.  In  the  Spectral  Data 
Transfer  System,  although  the  bounded  cross-correlation  characteristic 
is  not  being  utilized,  a  Gold  sequence  still  offers  more  security  than 
a  single  pseudorandom  sequence  because  there  are  more  codes  available 
for  altering. 

Two  pseudorandom-sequence  generators,  #1  and  #2,  of  the  same 
number  of  stages,  and  with  timing  pulses  provided  by  a  Bit-Rate  Clock, 
generate  two  sequences,  the  X  sequence  and  the  Y  sequence  respectively. 

The  X  sequence  and  the  Y  sequence  are  bit-by-bit  modulo-2  added  to  form 
a  Gold  sequence.  The  multiplexed  data  sequence  is  then  scrambled  by 
the  Gold  sequence  through  a  bit-by-bit  modulo-2  adder. 

The  frame-format  circuit  has  three  inputs:  the  scrambled  data 
sequence  and  the  X  sequence  both  at  the  bit  rate  and  the  frame  sync 
pattern  sequence  which  is  an  alternating  zero-one  sequence  at  ^  of  the 
bit  rate.  A  and  B  form  the  select  signal:  when  A=1  and  B=C,  the  output 
of  the  frame-format  circuit  is  the  framp  sync  L’t;  when  A=0  and  B=l, 
the  output  is  a  digit  from  the  X  seque  .ce;  and  at  the  remaining  K-2 
instances  both  A=0  and  B=0,  the  outputs  are  the  scrambled  data  bits.  Thus, 
the  output  sequence  of  the  frame-format  circuit  has  the  following  frame 
structure: 
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F1 

X. 

1 

S1 

pri 

Sk-2  | 

F2 

l+k  1 

LfL 

i**  i 

EJ 

Sk-2  ! 

F3 

.  .  . 

k  digits - -|- - -k  digits 


where  F 1  s  are  the  frame  sync  bits;  X.. ,  X.+^,...are  digits  from  the  X 
sequence;  and  s..  and  s.~*,  i  =  1,  2.....K-2  are  the  scrambled  data  bits. 

There  are  two  sequences  at  the  output  of  the  data  scrambler:  the 
scrambled  data  sequence  described  above  and  the  Y  sequence.  These  two 
sequences  are  to  be  interleaved  at  the  FEC  encoder.  Sync  recovery  will 
be  discussed  in  detail  in  later  sections;  it  suffices  to  indicate  here  that 
the  digits,  f°rm  a  sampled  X  sequence,  sampled 

every  K  digits,  are  used  to  recover  the  sync  of  the  X  sequence  at  the 
receiver  while  the  Y  sequence  transmitted  with  the  scrambled  data  is  for 
recovery  of  the  sync  of  the  corresponding  sequence. 
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3.3  Self-Synchronization 


At  the  receiving  end  the  descrambler  performs  the  reverse  process 
of  the  scrambler  as  shown  in  Figure  3.18.  Two  pseudorandom-sequence 
generators  identical  to  the  ones  in  the  scrambler  with  predetermined 
initial  states  are  required.  When  they  are  activated  at  the  right  moment, 
the  correct  Gold  sequence  is  generated  and  the  data  is  recovered.  However, 
in  case  synchronization  is  lost  means  must  be  provided  for  its  recovery. 

The  methods  of  sync  recovery  for  X  sequence  and  Y  sequence  are 
different.  As  mentioned  in  Section  3.2.4,  the  Y  sequence  is  being 
transmitted  along  with  the  scrambled  data.  The  received  Y  sequence  is 
stored  in  an  n-stage  buffer  register  and  is  compared  with  the  register 
contents  of  the  local  PN  generator  #2  through  a  bit-by-bit  comparator  as 
shown  in  Figure  3.3.  When  a  disagreement  is  detected  a  signal  is  sent 
to  replace  the  register  contents  of  the  local  PN  generator  #2  with  that  of 
the  received  Y  sequence.  The  disagreement  signal  also  activates  the  sync- 
recovery  process  for  the  X  sequence.  Thus  the  Y  sequence  performs  two 
functions:  sync  recovery  for  the  Y  sequence  and  sync  loss  detection. 

3.3.1  Preliminary  Description 

The  scrambled  data  sequence  is  the  bit-by-bit  modulo-2  sum  of 
the  data  sequence  and  the  Gold  sequence  which,  in  turn,  is  the  bit-by-bit 
modulo-2  sum  of  the  X  sequence  and  the  Y  sequence.  The  synchronization  of 
both  the  X  sequence  and  the  Y  sequence  is  required  for  descrambling  the 
data.  To  recover  the  sync  for  the  Y  sequence,  external  information,  besides 
the  scrambled  data  sequence,  is  needed;  in  our  case  the  external  information 
is  the  Y  sequence  itself  transmitted  along  with  the  scrambled  data  sequence. 
On  the  other  hand  no  external  information  is  needed  to  recover  the  sync 
of  the  X  sequence;  the  information  needed  is  embedded  in  the  scrambled  data 
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Figure  3.18  DESCRAMBLER 
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sequence.  This  is  called  "self-synchronization". 

Recall  the  frame  structure  of  the  scrambled  data  sequence. 

jh!  hi  b  I  ••  isi  •?  r,  j>lllti  q  i--|sf2  jV.’Wk  i  •••] 

J-*  -  k  digits - - -k  digits- - *-!-*■ - k  digits- - *- 

The  X's  in  the  scrambled  data  sequt-nce  are  from  the  original  X  sequence. 

pi-1  j  Xi.LXi+l  )  •  "  iXifk-2  )Xi  +  k-l  |Xi  +  k  jVk+1  j  ‘] Xi+1  k-1  ]Xi+2k|Xi+2k+l] 

—  -  k  digits-- . -  — . -k  digits - - — k  digits — 

At  the  receiver  these  digits 

V'-  x  i  +  k  *  *  *  x  i  +  2  k 

are  separated  from  the  scrambled  data  sequence.  When  sync  loss  is  detected 
n  of. these  digits 

xi *  xi tk  ’  xi +2k  *  ’ ' '  Xi  +  (n-l)k 

are  used  to  calculate  the  proper  initial  state  of  the  pseudorandom 
generator  for  sync  recovery.  The  logic  steps  of  the  sync  recovery  process 
is  shown  in  Figure  3.19. 

Before  we  attempt  to  explain  bow  this  self  synchronization  scheme 
works  we  shall  review  some  basics  of  Finite  Fields. 
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Figure  3.19  INITIAL  STATE  RECOVERY 


3.3.2  Pseudorandom  Sequence  and  Elements  of  a  Finite  Field 

It  has  been  described  in  Section  3.2.1  how  a  primitive  polynomial 
of  degree  n,  say,  P(x)  can  generate  a  pseudorandom  sequence  of  length 
L-2n-l.  The  same  primitive  polynomial  also  defines  a  finite  field  of 
2n  c-Ic'Vk  nts  and  we  want  to  describe  how  these  field  elements  are  defined. 

Consider  the  totality  of  polynomials  with  one  indetenninant  x  and 
with  coefficients  either  0  or  1.  Let  the  polynomials  be  represented  by 
F(x),  and  let  P(x)  be  a  primitive  polynomial  of  degree  n.  Divide  F(x) 
by  P ( x )  resulting  in 

F(x)  =  Q(x)P(x)  +  R(x) 

where  Q(x)  is  the  quotient  polynomial  and  R(x)  is  the  remainder  polynomial. 

The  remainder  R ( x )  is  of  degree  less  than  or  equal  to  r.-l; 

R(x)  r  cn_l  1  f  cn_2  x  2  +  •••  +  C}x  +  CQ 

where  is  either  0  or  1 .  There  are  altogether  2n  possible  remainder 
polynomials  and  they  form  a  finite  field  of  2n  elements,  sometimes 
referred  to  as  Galois  field  of  2n  elements  or  GF(2n)  for  short. 

Addition  and  multiplication  in  the  field  are  the  same  as  that 
hi  tween  .  . lynomials  except  that  in  the  case  of  multiplication  when  the  product 
is  of  degree  higher  than  n-1,  it  should  be  divided  by  P(x)  and  replaced 
by  the  remainder  polynomial.  For  example,  let 
F^x)  =  x4  r  x2  +  1, 

f  2 ( x )  =  +  *  +  1. 

and 

P(x)  x5  *  x2  +  1. 
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The  product  between  F^(x)  and  F2(x)  is 


F-j(x)  •  F^(x)  =  (x4  +  x2  +  l)(x3  +  x  +  1) 
7  4  2 

=  X  +  X  +  X  +  X  +  1  • 


Divide  F-j ( x )  •  F,,(x)  by  P(x)  resulting  in 


x7  +  x4  +  x2  +  x  +  1 


x2  (x5  +  x2  +  1)  +  (x  +  1) 


where  the  quotient  polynomial  is  x  and  the  remainder  polynomial  is 

c 

x  +  1.  Therefore  in  GF(2  )  the  product  between  the  element  F^x)  and 
the  element  F 2 ( x )  is  the  element  x  +  1: 

F-j  (x)  •  F2(x)  =  x  +  1. 

Since  a  polynomial  cun  be  put  into  a  one-to-one  correspondence 
with  its  coefficient  sequence,  i.e. 


F(x)  =  C  ,  x"'1  +  C  ,  xn'2  +...+C,x  +  C  *-*•  C 
n- I  n-2  I  0 


»  o>  •  •  >C-| ,  Cn, 


n-1 ’  n-2 


an  n-tuple  of  zeros  and  ones  can  be  regarded  as  an  element  in  GF(2n). 
Consider  the  example  in  Figure  3.9,  The  polynomial  P(x)  =x  +x  +lisa 
primitive  polynomial;  it  generates  a  maximum  length  sequence  of  length  L=31 
and  at  each  clock  time  it  generates  one  digit  of  the  sequence.  In  the 
mean  tine,  at  each  clock  time  the  contents  of  the  generator  are  a  5- 
tuple  which  can  represent  a  polynomial.  For  instance,  at  t=5,  the  output 
sequence  digit  is  "0"  and  the  field  element  is 
110  10 
which  represents  the  polynomial 
1  +  x  +  x3 
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and  at  t=24,  the  output-sequence  digit  is  "1"  and  the  field  element  is 

0  10  11 

which  represents  the  polynomial 

,3,4 
x  +  x  +  x  . 

There  are  31  distinct  5-tuples;  joining  all  the  zero  5-tuple,  00000, 
they  are  the  entire  32  elements  of  GF(2  ). 

We  have  established  a  one-to-one  correspondence  between  the 
digits  of  a  maximum-length  sequence  and  the  elements  of  the  field.  In 
the  example  shown  in  Figure  3.9  if  the  digit  of  the  output  sequence  at 
t ~0  is  identified  as  the  0-th  digit,  "0",  then  it  corresponds  to  the  field 
element  represented  by 
1  0  0  0  0 

and  the  5-th  digit  of  the  output  sequence,  "0",  corresponds  to  the 
field  element 

110  10 

and  the  24-th  digit  of  the  output  sequence,  "1",  corresponds  to  the  field 
element 

01011. 

Consider  again  the  example  shown  in  Figure  3.9.  The  polynomial 
generating  the  pseudorandom  sequence  is 
P(x)  =  x5  +  x2  +  1. 

The  reciprocal  polynomial  of  P(x)  is 
P*(x)  =  x5  +  x3  +  1. 

Connect  a  linear  feedback  shift  register  according  to  F*(x),  but  this 

time  the  modulo-2  adder  is  placed  between  the  shift  register  stages  as 

shown  in  Figure  3.20.  An  arbitrary  initial  state  of  a  5-tuple,  say  10011, 

is  loaded  into  the  shift  register;  after  30  shifts  the  register  contents 

returns  to  their  initial  state,  10011.  The  31  5-tuples  at  each  clock 
time  are  all  the  possible  5-tuples  except  the  all -zero  5-tuple, 

00000.  They  are  the  elements  of  GF(2  )  just  like  the  ones  in 
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Figure  3.9.  The  only  difference  between  these  two  sets  is  the 
order  in  which  the  elements  are  arranged. 

Let  us  make  two  important  observations. 

(1)  Compare  the  sequences  taken  from  the  5-th  stage  of  the  two 
shift  registers.  They  are  identical  except  for  the  point  at  which  each 
sequence  begins.  We  observe  that  these  two  sequences  are  identical  up 
to  a  phase  shift.  The  amount  of  phase  difference  depends  on  the  initial 
states  of  the  two  shift  registers. 

The  shift  register  in  Figure  3.9  is  sometimes  referred  to  as 
simple  shift-register  generator,  or  SSRG,  and  the  one  in  Figure  3.20  is 
referred  to  as  a  modular  shift-register  generator,  or  MSRG.  We  have 
illustrated  that  a  pseudorandom  sequence  generated  by  a  primitive 
polynomial,  P(x),  using  an  SSRG  is  identical  up  to  a  phase  shift  with 
that  generated  by  its  reciprocal  polynomial,  P*(x),  but  using  an  MSRG. 

(2)  Starting  from  the  5-tuple  1  0  0  0  0  in  the  list  of  the 
contents  of  the  shift  register  generated  by  P*(x)  =  x  +  x  +  1  using 
the  MSRG  in  Figure  3.20,  write  down  the  corresponding  polynomial  in 


ascending 

powers . 

1 

0 

0 

0 

0 

1  +  Ox  + 

Ox2 

+ 

Ox3 

+  Ox4 

0 

1 

0 

0 

0 

X 

0 

0 

1 

0 

0 

x2 

0 

0 

0 

1 

0 

X3 

0 

0 

0 

0 

1 

X4 

1 

0 

0 

1 

0 

1 

+ 

X3 

0 

1 

0 

0 

1 

X 

♦  X4 

1 

0 
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1 

0 

1  + 

X2 

+ 

X3 

0 

1 

0 

1 

1 

X 

+ 

3 

X 

+  x4 
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From  the  simple  calculations, 

X5  =  (1 ) ( 1  +  x3  +  X5)  +  (1  +  X3)  X5  =  1  +  X3 

x6  =  (x)  (1  +  x3  +  X5)  +  (x  +  X4)  X6  =  X  +  X4 

x7  =  (1  +  x2)(l  +  X3  +  x5)  +  (1  +  x2  +  x3)  X7  =  1  +  x2  + 

X8  =  (x  +  X3)(l  +  X3  +  x5)  +  (x  +  X3  +  X4)  «-#■  X8  =  X  +  X3  + 

it  can  be  verified  that  the  field  elements  are  arranged  in  the  order  of 
the  powers  of  >  ,  i .e. 

x1  (mod  P  (x)),  for  i  =  0,  1,  2,. ..,30. 


Now  i f  we 

let  B 

be  the  field 

element 

represented  by  the  polynomial 

x  or  the  n-tuple  0  1  0 

0  0,  then 

0 

1 

B  = 

B  = 

x 

2 

2 

B  = 

X 

„3 

3 

B  = 

X 

4 

4 

B  = 

X 

5 

1 

3 

=  1 

3 

B  = 

+  X 

+  B 

6 

4 

„  ^  4 

F  = 

X 

+ 

X 

B  +  B 

7 

1 

2  3 

=  1 

.  2  .  3 

B  = 

X  +  X 

+  +  B  +  B 

8 

3 

4 

„  .  „3  .  „4 

B  = 

X 

+  x  + 

x  = 

B  +  B  +  B 

Since  B8  =  1  +  B3, 

B  +  +  1  =  0 

i  .e.  P*(B)  =  B5  +  B3  +  1  =  0, 

which  means  B  is  a  root  of  P*(x).  Thus  we  have  reached  an  important 
result.  The  digits  of  a  pseudorandom  sequence  generated  by  a  primitive 


polynomial,  P(x),  of  degree  n  can  be  put  into  a  one-to-one  correspondence 

with  the  field  elements  of  GF(2n)  ordered  as  the  consecutive  powers  of 

a  root  6  of  P*(x),  which  is  the  reciprocal  polynomial  of  P(x).  Using 

5  2 

the  same  example  the  pseudorandom  sequence  generated  by  P(x)  =  x  +  x  +  1 

in  Figure  3.9  can  be  put  into  a  one-to-one  correspondence  with  powers 

5  3 

of  6  where  6  is  a  root  of  P*(x)  =  x  +  x  +  1.  This  correspondence  is 
shown  in  Figure  3.21. 


PSEUDORANDOM  SEQUENCE 
GENERATED  BY  P(x)  =  x5  +  x2  +  1 


0 

0 

0 

0 

1 

0 

1 

0 

1 

1 

1 

0 

1 

1 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

1 

1 

0 

1 

0 

0 
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POWERS  OF  6 

*.  .  c.  o 

p  (e)  =  r  +  +  l 

»°-i 

B1  =  8 
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B11 

B12 

B13 

B14 
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B16 

B17 

B18 

B19 

B20 

B21 

B22 

B23 

B24 

B25 

B26 

B27 

B28 

B29 

,30 


Figure  3.^1  CORRESPONDENCE  BETWEEN  A  PSEUDORANDOM 
SEQUENCE  AND  POWERS  OF  AN  ELEMENT  IN 
A  FINITE  FIELD 


3.3.3  Samp 1 ed_ Vers  ion  of  a  Pseudorandom  Sequence 

In  the  data  scrambler  system  described  in  Section  3.2,  two 
pseudorandom  sequences,  the  X  sequence  and  the  Y  sequence,  are  generated 
to  construct  the  Gold  sequence.  The  Y  sequence  is  transmitted  along  with 
the  scrambled  data  sequence.  In  every  k  digits  of  the  transmitted  data 
sequence  there  is  one  digit  which  comes  from  the  original  X  sequence. 

If  we  pick  out  this  digit  from  the  transmitted  sequence,  trie  resultant 
sequence  is  the  sampled  version  of  the  original  X  sequence,  sampled  every 
k  digits.  We  now  want  to  analyze  the  relationship  between  this  new  sequence 
and  the  original  X  sequence  and  to  find  out  how  this  new  sequence  can 
be  utilized  to  recover  the  synchronization  between  the  transmitted  X 
sequence  and  the  locally  generated  X  sequence  at  the  receiving  end.  For 
notational  convenience  in  the  following  discussion,  the  X  sequence  will 
be  called  the  "A  sequence"  and  the  new  sequence,  which  is  the  sampled 
version  of  the  A  sequence,  will  be  called  the  "B  sequence." 

Consider  the  pseudorandom  sequence  (A  sequence)  of  length  L=2n-1, 
which  is  generated  by  a  primitive  polynomial  of  degree  n,  P^(x). 

Sample  the  A  sequence  every  k  digits,  where  k  is  relatively  prime  to  L, 
and  form  the  B  sequence: 


k  +  1  ’ ' 


Since  k  and  L  are  relatively  prime  to  each  other,  the  B  sequence  will 
not  repeat  itself  until  every  digit  of  the  A  sequence  has  been  sampled; 
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hence  the  B  sequence  is  also  a  maxi  mum- length  sequence  and  the  length  of 
the  B  sequence  is  also  L. 

Since  the  B  sequence  is  a  maximum-length  sequence,  there  must 
be  a  primitive  polynomial,  Pg(x),  which  will  generate  the  B  sequence. 

By  the  important  result  obtained  in  Section  3.3.2,  that  the  digits  of 
a  pseudorandom  sequence  generated  by  a  primitive  polynomial,  P(x),  of 
degree  n  can  be  put  into  a  one-to-one  correspondence  with  the  field 
elements  of  GF(2n)  ordered  as  the  consecutive  powers  of  a  root  of  the 
reciprocal  polynomial,  P*(x),  of  P(x)  the  A  sequence  corresponds  to 
powers  of  a  *  and  the  B  sequence  corresponds  to  powers  of  b~*: 


where  a'1  is  a  root  of  Pj(x)  and  p"1  is  a  root  of  P*(x)  and  p'1  =  (a'1)k. 

In  Section  3.2.1  we  have  defined  the  reciprocal  polynomial  without 
referring  to  field  elements  as  follows:  if  P ( x )  is  a  polynomial  of  degree 
n,  then  the  reciprocal  polynomial 
P*(x)  *  xnP(*-). 

In  terms  of  field  elements,  reciprocal  polynomial  has  the  following  meaning 
if  a  field  element,  Y,  is  a  root  of  a  polynomial  P(x),  i.e. 


is  a  root  of  the  reciprocal  polynomial 


p(y)  =  o 

then  the  inverse  element  of  Y,  Y'^, 

P*(x)  or 

P*(y-1)  =  0. 

5 

As  an  example,  consider  the  finite  field  GF(2  )  generated  by 
primitive  polynomial 

P(x)  =  x5  +  x3  +  1 
the  field  element 


Y  =  0  1  0  0  0 

is  a  root  of  P(x),  i.e.  P(y)  =  0. 

This  can  be  verified  as  follows:  from  Figure  3.20,  we  observe  that 

Y5  =  (1  0  0  1  0) 

Y3  =  (0  0  0  1  0) 

_ 1  =(10000) 

Y5  +  y3  +  1  =  (0  0  0  0  0) 

Since  the  5-tuple,  0  0  0  0  0,  is  the  field  element,  0,  therefore 
P(y)  =  y5  +  y3  +  1  =  0. 

The  reciprocal  polynomial  P*(x)  is 
P*(x)  =  x5  P(£) 


5  2 

=  x3  +  x^  +  1 


and 
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(Y-1)5  =  Y"5  =  >26  =  (1  1  0  1  0) 

(y-1)2  =  r'2  =  y29  =  (o  i  o  i  o) 

1  _ _ _ (1  0  0  0  0) 

(y-1)5  +  (Y-1)2  +1  =(00000) 

hence 

P*(y_1)  =  (y'1)5  +  (y'1)2  +1=0 

Now  from  the  results  obtained  above, 

-1  ,  -1 ,k 

6  -  (a  )  , 

since  each  field  element  has  a  unique  inverse,  it  follows  that 


Furthermore  a-1  is  a  root  of  P*(x)  and  e"1  is  a  root  of  Pg(x)  imply 
that  a  is  a  root  of  P^( x)  and  6  is  a  root  of  Pg(x)  respectively.  It 
can,  therefore,  be  concluded  that  if  a  is  a  root  of  a  primitive 

polynomial,  PA(x),  generating  the  A  sequence  and  g  is  a  root  of  Pg(x) 

generating  the  B  sequence;  and  if  the  B  sequence  is  a  sampled  version 

of  the  A  sequence,  sampled  every  k  digits  then  B  is  the  k-th  power  of 

a,  i .e. 
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When  k  =  2m,  for  m  £  n-1,  pg(x)  =  P^(x).  This  follows  from  the 
fact  that  in  GF(2n)  if  a  is  a  root  of  a  polynomial  P(x),  then  a2,  aA,..., 
2^ 


a  are  also  roots  of  P(x).  The  notation  Pg(x)  =  P^(x)  means  that  the  B 
sequence  is  identical  to  the  A  sequence  up  to  a  phase  shift  [7]. 
Consider  the  following  examples. 


Example  1 . 


n  =  3,  L  =  2-1  =  7 
.3 


Choose  Pft(x)  =  x  +  x  +  1  with  1  0  1  as  the  initial  state  as  in  Figure 
3.7.  The  generated  A  sequence  is: 


\ 

Ao 

A1 

A2  A3 

A4 

A5 

A6 

1 

0 

1  0 

0 

1 

1  . 

t 

If  k  =  22 

=  4, 

then 

i 

B0 

=  Ao 

=  1 

B4  =  A16=  A2  =  1 

B1 

II 

4^ 

=  0 

B5  =  A20=  A6=  1 

i 

{ 

I. 

B2 

A8 

=  A1  = 

0 

B6  =  A24  =  A3  =  0 

1; 

l 

B3 

=  A12  A5 

=  1. 

The  B  sequence 

is 

;{ 

B0 

B1 

B2 

B3 

B4 

B5 

B6 

1  i 

1  t 

1 

0 

0 

1 

1 

1 

0  . 

Comparison  between  these  two  sequence  shows  that  the  sampled  sequence  is 
same  as  the  original  sequence  except  for  the  starting  digits. 

If  k  =  3,  then 


B0 


1 


B4  “  A12  =  A5  "  1 


Bi  =  a3  -  0 


B5  "  A15  =  A1  "  0 


Bo  =  Ae 


«  A18=  A,  =  0 


B3  =  Ag  -  A2  -  1. 
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The  B  sequence  is: 


B2  B3  B4 


This  sequence  is  different  from  the  original  A  sequence.  Since  there 
are  only  two  different  pseudorandom  sequences  for  n=3,  the  polynomial 
PB(x)  which  generates  the  B  sequence  must  be  the  reciprocal  polynomial 
of  PA(x);  hence 

PB^  =  PA^  x  +  x  +  1  • 

Example  2. 

n  =  4,  L  =  24-  1  =  15 

Choose  P^U)  =  x4  +  x  +  1  with  1  1  0  1  as  the  initial  state  as  in 
Figure  3.8.  The  A  sequence  is 


A0 

A1 

A2  A3 

A4 

A5 

A6 

A7 

A8 

A9 

A10 

A11 

A1 2 

A1 3 

A14 

1 

0 

1  1 

0 

0 

1 

0 

0 

0 

1 

1 

1 

1 

0  . 

For  k 

-  22 

=  4, 

,  the  B 

sequence 

is 

B0 

B1 

82  B3 

B4 

B5 

B6 

B7 

B8 

B9 

B10 

B11 

B12 

B1 3 

B14 

1 

0 

0  1 

0 

0 

0 

1 

1 

1 

1 

0 

1 

0 

1 

which 

is  a 

shifted  version 

of 

the 

A  sequence. 

For  k 

=  7, 

the 

B  sequence 

is 

B0 

B1 

B2  B3 

B4 

B5 

B6 

B7 

B8 

B9 

B10 

B11 

B1 2 

B1 3 

B14 

1 

0 

0  1 

1 

0 

1 

0 

1 

1 

1 

1 

0 

0 

0 

This  sequence  is  different  from  the  original  A  sequence.  Since  there 
are  only  two  distinct  pseudorandom  sequences  for  n  =  4,  the  polynomial 
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which  generates  the  B  sequence  must  be  j 

PB(x)  *  PA(x>  "  x4  +  x3  ♦  1.  ! 

Example  3. 

n  =  5,  L  =  25  -  1  =  31 

Choose  P ^(x)  =  x^  +  x^  +  1  with  1  0  0  0  0  as  the  initial  state  as  in 

Figure  3.9.  The  A  sequence,  the  B  sequence  by  sampling  every  4-th  digit, 

and  the  B  sequence  by  sampling  every  5-th  digit,  are  shown  in  Figure  3.22. 

2 

Again,  the  B  sequence  for  k  =  2  =  4  is  identical  with  the  A  sequence 

up  to  a  phase  shift  and  the  B  sequence  for  k  =  5  is  different  from  the  i 

A  sequence.  From  the  Table  of  Irreducible  Polynomials  we  found  the 

! 

polynomial 

j 

5  67H- 

Thus ,  for  k  =  5  j 

PB{x)  =  x5  +  x4  +  x2  +  x  +  1.  | 

1 

i 

j 

1 

We  now  summarize  the  results  as  follows.  Let  Pfl(x)  be  a  primitive 
polynomial  of  degree  n  which  generates  a  pseudorandom  A  sequence  of 
length  L  =  2n  -  1.  Sample  the  A  sequence  at  every  k-th  digit;  when 
(K,L)=1,  the  sampled  sequence  (the  B  sequence)  is  also  pseudorandom  and 
is  generated  by  a  primitive  polynomial  P&(x) .  If  k=2m,  for  m  <  n-1,  then 

pB(x)  =  PA(*) 

which  means  that  the  B  sequence  and  the  A  sequence  are  identical  up  to  a 
phase  shift.  If  k  1  2m,  then 

PB(X)  t  PA (*)• 
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A  SEQUENCE 

B  SEQUENCE 
k  =  4 

B  SEQUENCE 
k  =  5 

A0 

0 

B0 

0 

B0 

0 

A1 

0 

B1 

1 

B1 

0 

A2 

0 

B2 

1 

B2 

1 

A3 

0 

B3 

1 

B3 

0 

A4 

1 

B4 

0 

B4 

1 

A5 

0 

B5 

1 

B5 

1 

A6 

1 

B6 

1 

B6 

1 

A7 

0 

B7 

0 

B7 

1 

A8 

1 

B8 

0 

B8 

1 

A9 

1 

B9 

0 

B9 

0 

A10 

1 

B10 

1 

B10 

1 

A11 

0 

B11 

1 

B11 

1 

A12 

1 

B12 

1 

B12 

0 

A13 

1 

B13 

1 

B13 

0 

A14 

0 

B14 

1 

B14 

1 

A15 

0 

B15 

0 

B15 

1 

A16 

0 

B16 

0 

B16 

1 

A17 

1 

B17 

1 

B17 

0 

A18 

1 

B18 

1 

B18 

0 

A19 

1 

B19 

0 

B19 

0 

A20 

1 

B20 

1 

B20 

0 

A21 

1 

B21 

0 

B21 

1 

A22 

0 

B22 

0 

B22 

1 

A23 

0 

B23 

1 

B23 

0 

A24 

1 

B24 

0 

B24 

1 

A25 

1 

B25 

0 

B25 

0 

A26 

0 

B26 

0 

B26 

1 

A27 

1 

B27 

0 

B27 

0 

A28 

0 

B28 

1 

B28 

0 

A29 

0 

B29 

0 

B29 

1 

o 

co 

<r 

1 

B30 

1 

B30 

0 

Figure  3. .22  PSEUDORANDOM  SEQUENCE  AND  ITS  SAMPLED  SEQUENCES 
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However,  if  c 

is  a  root  of 
of  the  Table 


is  a  root  of  P^(x),  then 

a  k 
e  =  a 

PB(x).  The  polynomial  Pg(x)  can  be  determined  with  aid 
of  Irreducible  Polynomials. 


90 


V\ 


I 


» 


t 


I 


I 


3.3.4  Central  Idea  of  the  Self- Synchronization  Scheme 

The  A  sequence,  which  is  the  original  X  sequence,  is  a  component 
of  the  Gold  sequence.  At  the  receiving  end  an  identical  pseudorandom- 
sequence  generator  generates  the  A  sequence  in  order  to  perform  the 
descrambling  operation.  When  a  sync  loss  is  detected,  the  contents  of  that 
generator  must  be  cleared  and  the  "correct  initial  state"  must  be  loaded 
into  the  shift  register  to  recover  the  synchronization.  The  problem  is 
how  to  determine  the  correct  initial  state. 

Every  k-th  digit  of  the  A  sequence  is  transmitted  along  with 
the  scrambled  data  sequence  and  is  received  at  the  receiving  end.  The 
sequence  of  these  received  digits  is  called  the  B  sequence  which  is  the 
sampled  version  of  the  A  sequence.  Suppose  that  at  the  time  sync  loss 
is  detected  one  such  B  digit  is  just  received.  We  label  this  digit  as 
Bq,  corresponding  to  A g  in  the  A  sequence.  The  (n-1)  succeeding 
digits,  A^  A2,...,Anl,  together  with  AQ  are  the  correct  initial  state 
beginning  from  Ag.  These  digits,  Aj,  can  be  obtained  from 

the  B  sequence.  Specifically,  since  the  sample  interval,  k,  and  the 
length  of  the  A  sequence,  L,  are  relatively  prime  to  each  other,  the 
B  sequence  will  not  repeat  itself  until  every  digit  in  the  A  sequence  is 
sampled  and  this  occurs  after  exactly  k  cycles  of  the  A  sequence.  For 
instance,  if  L  =  7  and  k  =  3  the  B  sequence  will  repeat  after  3  cycles  of 
the  A  sequence. 


- — Hh 

Ac  A- 


A0  A1  A2  A3  A4  A5  A6  Ap  A1  A2  A3  A4  A5  A6  A0  A1  A2  A3  A4  A5  "6  "0 


Now,  there  must  be  a  B-digit  which  equals  4-j .  Let  this  B-digit  be 
kp  =  1  {mod  L). 


Bp;  then 
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This  is  equivalent  to  measuring  the  same  distance  with  two  different 
yardsticks:  p  multiples  of  k  must  be  equal  to  m  multiples  of  L  plus  one. 
Therefore, 


A0  B0 


Aj  =  Bp  ,  where  kp  =  1  (mod  L), 


It  is  not  difficult  to  show  that 


2p 


A  B/  . ,  . 
n-1  ( n- 1 } p 


However,  this  is  not  acceptable  because  it  takes  too  long  to 
wait  for  all  these  digits.  To  have  some  idea  as  to  how  long  it  will 
take  to  accumulate  these  digits,  let  us  consider  an  example.  Suppose 
the  PN  generator  has  11  stages,  i.e 


n  =  11. 

Then 

L  =  211  -  1  =  2047. 

Let 

k  =  13. 

Solve 

kp  =  1  (mod  2047) 

giving 

p  =  315, 

for  13 

x  315  =  4095  =  2  x  2047 

Therefore 

A0  =  B0 


A1  B  31 5 


A2  "  B630 


A3  '  B945 


A4  "  8 1 260 


A5  =  B 1 575 


A6  =  8 1390 
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A7  =  B2205  =  B 1 58 
A8  =  B2520  =  B473 
A9  =  B2335  =  B  788 
A10  =  B3150  =  B 1 103  ‘ 


We  need  to  wait  for  1890  B-digits  which  is  equivalent  to  1890  x  13  =  24,570 
A-digits.  This  means  that  when  we  lose  sync,  24,570  digits  later  we 
will  find  the  current  initial  state  beginning  at  Ag.  If  we  want  to 
recover  the  sync,  we  have  to  speed  up  the  shift  register  to  catch  up 
these  24,570  digits. 

The  central  idea  of  the  self  synchronization  scheme  is  to  reduce 
this  waiting  time  to  less  than  a  bit  duration.  When  the  Spectral  Data 
Transfer  System  is  running  in  a  normal  mode,  the  scrambled  data  sequence 
is  being  received  and  being  descrambled;  in  every  frame  of  k  bits  one 
B-digit  is  also  received.  The  B-digits  are  loaded  serially  into  an  n-stage 
shift  register;  thus  n  B-digits  are  being  stored  and  up-dated  at  all  times - 
When  sync  loss  is  detected  the  most  recent  digit  entering  the  shift 
register  is  designated  as  Bn  and  the  one  ready  to  be  discarded,  as  BQ. 
The  digit  BQ  =  Ag  occurred  (n-l)k  digits  ago  in  the  A  sequence.  The 
self-synchronization  scheme  is  to  determine  the  initial  state,  Ag,  A^, 
A2,...,An  j,  beginning  at  Ag  from  the  presently  available  B-digits, 

Bq,  Bj,  B2,...,Bn_r 


Take  the  same  example  of  n  =  11. 
Sync  Loss 
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It  can  be  seen  from  this  example  that  if  we  can  determine  the 
initial  state  A0,  Aj,...,Ajq  beginning  at  Aq  immediately  at  the  instant 
sync  loss  is  detected  we  have  130  digits  to  catch  up  instead  of 
24,570  digits. 

Since  it  has  been  recognized  that 


A0  B0 


Al  =  Bp 


A2  “  B2p 


An-1  B(n-l)p 


the  problem  is  now  to  calculate  Bp,  B2p,  .  ...B^  ^p  in  terms  of  BQ,  Bj 


B2* ‘ ' ,Bn-l‘ 


n-1 


3.3.5  Determination  of  Bm  =  b-B-  for  all  m  <  L 

m  i =0  1  1 


It  has  been  establ ished  in  Section  3.3.2  that  the  digits  of  a 
Dseudorandom  sequence  generated  by  a  primitive  polynomial  of  degree  n, 
say  Pg(x),  can  be  put  into  a  one-to-one  correspondence  with  the  powers 
of  an  element  of  GF(2n)  and  that  this  field  element  is  a  root  of  the 


reciprocal  polynomial  of  Pg(x)  denoted  by  PR(x).  For  example,  if 


B' 


5  2  5  2 

PB(x)  =  x  +  x  +  1  the  pseudorandom  sequence  generated  by  x  +  x  +  1 


is  shown  in  Figure  3.9.  The  reciprocal  polynomial  Pg(x)  is  x5  +  xJ  +  1 


Powers  of  p,  Pg(B)  =  0,  are  generated  in  Figure  3.20.  Since 


{$  =  1  0  0  0  0 


6  =  01000 

B2  =  0  0  1  0  0 

63  =  0  0  0  1  0 

e4  =  o  o  o  o  i 
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and  any  power  of  g  is  an  n-tuple,  gm  can  be  expressed  as  a  linear 
0  2  3  4 

combination  of  g  ,  g,  g  ,  g  ,  and  g  ,  i.e. 


m 


£ 

T^O 


b,  g' 


On  the  other  hand,  Bm  is  just  the  n-th  component  of  the  n-tuple  representing 
gm;  therefore 

n-1 

m  \  '  ,  i 
6  ■  bi  6 
i=0 


implies  [7] 


n-1 


i=0 


b. 

l 


B.. 

i 


For  example  in  Figure  3.20 
g19  =  1  1  0  1  1 

=  10000  +  01000  +  00010  +  00001 
=  1  +  6  +  g3  +  g4 


and  in  Figure  3.9 


B0  -  0 

Bj  =  0 

b3  -  0 

b4  -  1 

fi19  =  B0  +  B1  +  B3  +  B4 


=  0  +0  +0  + 


1. 


Moreover,  in  the  pseudorandom-sequence  generator  using  the  simple 
shift-register  generator  the  consecutive  digits  of  the  sequence  shift 
one  stage  down  from  left  to  right  at  each  clock  pulse;  hence  Bm  can  be 
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obtained  with  a  modulo-2  adder  and  proper  connections  as  shown  in 
Figure  3.23.  At  the  next  clock  pulse  6^  replaces  B4,  B4  replaces  63,..., 


Sj  replaces  Bq,  and  the  output  of  the  modulo-2  adder  is 
B5  +  B4  +  B2  +  B1  B20' 

The  modulo  2  adder  and  the  connections  are  sometime  called  the  phase 
shift  network. 

It  should  be  noted  that  this  property  is  not  applicable  to  the 
modular  shift-register  generator  with  the  reciprocal  polynomial  although 
the  output  at  the  last  stage  produces  the  identical  pseudorandom 
sequence. 

We  are  now  able  to  determine  A,  which  is  equal  to  B  , 
n-1  1  p 


B 


■£bi  v 


i=0 


Let  us  consider  the  following  examples. 
Example  Is.. 

n  =  3 

PA(x)  =  x3  +  x  +  1 
(i)  k  =  22  =  4 

pB(x)  =  PAU)  =  x3  +  x  +  1 
and  p  =  2  =2 


The  circuit  to  calculate  Bp  is  shown  in  Figure  3.24. 


(ii)  k  =  3 

PB(X)  =  X3  +  X2  +  1 
Solve  for  p  in 
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kp  =  1  (mod  7) 
p  =  5 

P*(x)  =  x3  +  x  +  1 
Pg( b)  =  b3  +  b  +  i  =  o 
B3  =  B  +  1 

B4  =  B  B3  =  B(B  +  1)  =  B2  +  B 

S  4  ?  o  o 

B  =  B  B  =  B(6  +  B)  =  B  +  B 
=  B  +  1  +  B2 
=  B2  +  B  +  1 
B5  =  B2  +  B1  +  B0 

The  circuit  for  calculating  Bp  is  shown  in  Figure  3.25. 
Example  2. 

n  =  4 

PA(x)  =  x4  +  x  +  1 

(i)  k  =  22  =  4 

Pg(x)  =  PA(x)  =  X4  +  X  +  1 

p  =  2n-m  =  4 

*,  ,  4  3 

PB(X)  =  X4  ♦  xJ  ♦  1 

*.  .  4  3 

Pg(  B)  =  B  +  B  +  1  =  0 

64  -  S3  ♦  1 
Bp  ‘  S4  "  B3  +  B0 

The  circuit  for  calculating  Bp  is  shown  in  Figure  3.26. 

(iii)  k  =  7 

PB{X)  .  x4  ♦  x3  *  1 

7p  =  1  (mod  15) 

p  =  13  [13  x  7  =  91  =  6  x  15  +  1] 


_ >  v N 


★  ,  4 

PB(x)  =  X  +  X  +  1 
Pr(b)  *  B4  +  B  +  1  =  0 


B  =  B  +  1 

B5  =  B2  +  B 
B6  B3+B2 
B7  =  B4  +  B3 
=  B3  +  B  +  1 
B8  =  B4  +  B2  +  B 


=  8+1 
B9  =  B3  +  B 

B10  =  B4  +  B2 


=  6  +  B  +  1 

11  3  2 

B  =  B  +  B  +  B 

12  4  3  2 

B  =  B  +  B  +  B 

3  2 

=  B  +  B  +  B  +  1 

13  4  3  2 

B  =  B  +  B  +  B  +B 

3  2 

=  BJ  +  B  +  1 

B  =  B.,  =  B-  +  B9  +  Bn 
p  13  3  2  0 

The  circuit  for  calculating  Bp  is  shown  in  Figure  3.27. 


Example  3. 


n  =  5 

Pfi(x)  =  x5  +  x2  +  1 
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(i)  k  =  2  =  4 


PB(x)  =  PA(x)  =  x5  +  x2  +  1 


p  =  2n'm  =  25-2  =  8 


Pg(x)  =  X5  +  X-5  +  1 


p*(b)  -  b5  +  b3  +  l  »  o 


b5  -  b3  «■  i 


b6  =  b4  +  e 


b7  -  b5  +  b2 


=  63  +  32  +  1 


b8  =  b4  +  b3  +  b 


8R  =  B4  +  +  B1 


The  circuit  for  calculating  Bp  is  shown  in  Figure  3.28. 


(ii)  k  -  5 


PB(X)  *  X5  ♦  X4  ♦  X2  +  X  +  1 


5o  =  1  (mod  31) 


p  =  25  [5  x  25  =  125  =  4  x  31  +  1] 

* 

V 

* 

V 


*  B  4  3 

PR(x)  =X  +  X  +X  +X+1 


P*(B)  =  b5  +  B4  +  B3  +  B  +  1  =  o 


c  a  3 
B  =  B  +  B  +  B  +  1 


B6  =  B5  +  B4  +  B2  +  B 


3  2 

=  8  +  B  +  1 


B12  =  B6  +  B4  +  1 


4  3  2 

B  +  8  +  8 


24  .  e8  *  e6  *  e 4 
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=  5,  k  =  4) 


_ ^s.-  .**  V  \ 


=  b2(b3  +  b2  +  l)  +  b3  +  b2  +  i  +  b4 

4  3  4  2  3  2  4 

=  B  +  B  +  B  +  1  +  B  +  B  +  B  +  B  +  1  +  B 

4 

=  B  +  B 

25  4  3  2 

B  =  B  +  B  +  B  +  1  +  B 

4  3  2 

=  B  +  B  +  B  +  B  +  1 

Bp  -  B25  “VVVV'o 
The  circuit  for  calculating  is  shown  in  Figure  3.2?. 

3.3.6  Calculation  of  B^p,  B^,. ..  ,B^n_^p 

In  Section  3.3.5  we  have  determined  the  coefficients  bi  for 
calculating 

n-1 

B  =  /  _,  b.B.  for  all  m  <  L. 

m  x  .j  li  — 

i=0 

and  we  have,  in  particular,  found  the  coefficients  for  Bp.  These 

coefficients  enable  us  to  construct  a  connection  network,  which  is 
called  the  phase-shift  network,  to  obtain 


directly  from  the  sampled  digits  of  the  original  pseudorandom  sequence. 

In  principle  we  can  use  the  same  procedure  to  construct  a  phase 
shift  network  for  each  component  of  the  initial  state.  However,  a  close 
observation  reveals  that  the  connection  network  already  obtained  can 
calculate  the  remaining  components,  -  ^2p’‘''’^n-l  =  ^(n-l)o’  as  we^ * 
G1ven  n-l 

Bp  *  £  biBi 

i=0 
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which  corresponds  to 
n-1 


■  E  bi  * 


i=0 


where  Bp  is  the  p-th  digit  of  the  pseudorandorr.-sequence  generated  by 

the  polynomial  Pg(x)  and  BP  is  the  p-th  power  of  a  field  element,  B, 

which  is  a  root  of  P£(x),  the  reciprocal  polynomial  of  P  (x),  B„  can 

**  B  2p 

be  calculated  as  follows: 

e2p  =  6p  •  eD  =  6p  bt  8' 


i=0 


n-1 


E  bi^p  +  1 


i=0 


which  corresponds  to 
n-1 

B2p  ’  E  »i 

i=0 


P+i  ' 


Using  a  similar  derivation,  we  obtain 
n-1 


B 


3.  =  ^  b. 
JP  Z-/  i 


i  =  l 


B(j-l)p+i,  j 


=  2,3, _ ,n-l. 


We  now  want  to  show  how  B.  can  be  calculated  through  a  circuit 

J  r 

implementation.  Let  us  consider  the  following  examples. 

Example  1. 

n  =  3,  k  =3,  p  =  5 
P8(x)  »  x3  ♦  x2  +  1 

and  Bp  =  B5  ■  B2  +  Bj  +  B„ 

Construct  the  pseudorandom-sequence  generator  and  the  phase  shift  network 
as  in  Figure  3.30*  The  digits  Bg,  ,  and  are  in  the  shift  register  at 


'  t 
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the  instant  when  sync  loss  is  detected.  Then 


VB0 


and 


=  Bg  throught  the  connection  network. 
fl2  ‘  B2p  =  B10  "  B5+2  +  B5+l  +  B5 


“  B7  +  B6  +  BS 


These  digits  B^,  B^,  and  B,-  are  obtained  at  the  output  of  the  phase- 
shift  network  after  each  shift  of  the  pseudorandom  sequence  generator. 
Example  2 


For  n  =  4  and  k  =  7. 


Pg(x) 


4  3 

x4  +  x  +  1 


kp  =  1  (mod  2n-l)  =t>  p=13 
and  Bp  =  Bl3  =  B3  +  B2  +  Bq 

B2p  =  B26  =  Bp+3  +  Bp+2  +  BP  =  B16  +  B15  +  B13 

or  B2p  =  B26  =  B11  =  B1  +  B0  +  B13 


Similarly, 


B3p  B39  B9  "  B14  +  B13  +  B 


11 


The  principle  of  operation  is  illustrated  in  Figure  3.31. 
Example  3 


For  n  =  5  and  k  =  6 
PB(x)  =  x5  +  x4  +  x3  +  x  +  1 
p  =  25 

and  Bp  =  B^  =  B^  +  B3  +  B^  +  B^  +  Bg 

B2p  =  B50  =  B19  =  B29  +  B28  +  B27  +  B26  +  B25 

B3p  =  B75  =  B13  =  B23  +  B22  +  B21  +  B20  +  B19 

B4p  =  B100  =  B7  =  B17  +  B16  +  B15  +  B14  +  B13 

The  principle  of  operation  is  illustrated  in  Figure  3.32. 
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Following  the  above  examples,  we  can  describe  the  principle 
of  the  operation  of  the  essential  part  of  the  sync  recovery  circuit  as 


follows . 

The  receiver  receives  and  stores  n  sampled  digits  of  the  original 

pseudorandom  sequence  in  an  n-stage  shift  register;  at  the  instant  when 

sync  loss  is  detected,  the  feedback  connections  corresponding  to  Pg(x) 

are  activated  and  the  shift  register  becomes  a  pseduorandom-sequence 

generator  generating  the  B  sequence  which  is  the  sampled  sequence 

corresponding  to  original  sequence.  At  the  same  time  the  phase  shift 

network  is  also  activated.  The  values  of  Art=Bri'and  A.=B  are  thus  ob- 

u  u  l  p 

tained  immediately. 

At  each  clock  time  another  dioit,  B  ...  appears  at  the  output 

D+l 

of  the  phase-shift  network.  These  diqits.  B  ,  B  ...... B  .  .  are  to 

p  p+1  p+n-1 

be  stored  in  an  n-stage  shift  register  buffer.  At  the  (n-l)-st  clock 
time  the  buffer  is  filled  up,  the  contents  of  the  pseudorandom  sequence 
generator  Pg(x)  are  cleared,  and  replaced  by  the  contents  of  the  buffer. 
The  value  of  A2=^2p  now  aPPears  at  the  output  ot  the  phase  shift 
network.  Repeat  this  procedure  n-2  times  in  order  to  obtain  the  entire 
initial  state: 

Mq,  Aj,  A2,. . . ,An_j. 


Load  this  initial  state  into  the  pseudorandom  sequence  generator  P.(x) 

rt 

of  the  original  pseudorandom  sequence,  and  speed  up  the  clock  to 
recover  the  sync. 


3.3.7  Determination  of  the  Phase  Shift  Network 

The  connections  of  the  phase  shift  network  are  constructed 
according  to  the  coefficients  b..  in 
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■  Z 


Since  8  is  a  root  of  the  reciprocal  polynomial  Pg(x), 
pg(e)  =  0 

and  the  degree  of  the  polynomial  Pg(x)  is  n  and  8n  is  a  linear  combination 
of  the  terms  of  degree  less  than  or  equal  to  (n-1).  Therfore  6^  can  be 
calculated  by  successive  application  of  the  recursion  formula 
SJ+1  =  e-8J 

and  successive  reduction  of  the  degree  of  the  term  8n. 

This  procedure  is  simple  in  concept,  but  very  lengthy  and  tedious 
in  practice,  and  almost  impossible  if  n  is  very  large.  We  shall 
present  a  much  shorter  method  by  considering  the  following  two  examples. 

Example  1. 

n  =  11,  k  =13 

Choose  P^U)  =  x11  +  x^  +  1. 

Since  k  =  13  /  2m,  for  m  <  n  =  11 
PB(x)  f  PA(x) 

and  Pg(x)  has  a  root  which  is  the  13-th  power  of  a  root  of  pA(x).  From  the  Table 
of  Irreducible  Polynomials,  it  is  found  that 
pB(x)  =  x11  +  x6  +  x5  +  x  +  1  . 

The  reciprocal  polynomial  of  Pg(x)  is 

P*(x)  =  x11  +  x10  +  x6  +  x5  +  1. 


Next,  solve  for  p  from 


kp  =  1  (mod  2047),  L  =  2l  -  1  =  2047 
13  x  315  =  4095  -  2  x  2047  +  1 


yielding  p  =  315. 


Hence 


A0  =  B0 
A1  =  B315 ' 


We  need  to  find  an  expression  for  in  terms  of  Bq,  Bj,  B^.-.-jB^g. 


>  *  S 


The  next  step  is  to  write  a  computer  program  simulating  the 
linear  feedback  shift  register  in  an  MSRG  configuration  according  to 

Pg(x)  =  x11  +  x10  +  x3  +  x3  +  1  as  shown  in  Figure  3.33  with 

b°  =  10000000000 

5  =  01000000000 
B2  =00100000000 

510  00000000001 

315 

and  calculate  8  according  to 


(a) 

(b) 


Calculate  B11,  B12,...,B3y  resulting  in 
B39  =  bQ  +  bjB  +  b2B2  +...+  b10810. 

Calculate  839,  B39+1,  B39+2 . B39+1°  and 


„78  =  k  q39  +  u  Q 39+1  .  .  39+2  . 

8  Uq8  +  bjB  + 


39+10 

8 


(c)  Repeat  (b)  until 

Q 312  _  ,  4x39  ,  .  4x39+1  ,  4x39+2, 

8  =  Dg6  +  OjB  +  D^B  +...+D 

(d)  Calculate  B3*3,  B3^,  B3*3  obtaining 

B315  =  b10  +  B7  +  B2  +  B  +  1- 


The  sync  recovery  circuit  is  constructed  as  in  Figure  3.34. 


4x39+10 

10P 


Example  2. 

n  =  39,  k  =  8 

Choose  P^(x)  =  x39  +  x^  +  1 
Since  k  =  8  =  23 
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Figure  3.33  MODULAR  SHIFT  REGISTER  GENERATOR, 


Figure  3.34  SYNC  RECOVERY  CIRCUIT 


PB(x)  =  PA(x)  =  X39  +  X4  +1. 

The  reciprocal  polynomial  of  Pg(x)  is 

Pg(x)  =  x39  +  x35  +  1. 

Solve  for  p  from 

kp  =  1  (mod  239  -  1) 

giving  p  =  239-3  =  238 

23^  0  12  38 

We  need  to  determined  in  terms  of  8  ,  8  ,  B  6  .  First  we 

write  a  computer  program  simulating  the  modular  shift  register  for 

P*(x)  as  shown  in  Figure  3.35  with 

8°  =  1  0  0  0 . 0 

6  =  010  0 . 0 

62  =  0  0  1  0 . 0 

638  =  0  0  0  0 . 01. 

' - ' 

39  digits 

The  calculation  procedure  is  as  follows. 

(a)  Calculate  63^,  84°,...,B64  and 

6^4  =  1>q  +...+t>2g33^' 

(b)  Calculate  B64,  B64+1,B64+2,. . . ,B64+3S  and 

B128  =  B2  =  b/4  ♦  b/4*1  ♦  f>/4+2  t...t638B64+3S. 

(c)  Repeat  (b)  . 

/.l/n,/*1*!/"  +...+b38B27+38 

„36  ?35  ?35+i  735  „  935 

62  =  b0BZ  +  b^2  1  +  b2e2  +4+...+b38B2  +38 

(d)  The  solution  is 

238  24  22  5 

6  =  8  +  8  +  e  . 
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Therefore  the  sync  recovery  circuit  is  constructed  according  to  Figure 
3.36.  From  the  two  examples  we  can  appreciate  the  powerfulness  of  the 
method,  especially  in  Example  2.  We  can  roughly  estimate  the  number  of 
calculations: 

initial  calculation  (step  a)  =  64  -  38  =  26 

+  (36-6)  x  38  =  1,140 

Total  1,166. 

Otherwise  the  number  will  be  approximately  2 ^  =  6.8  x  101Q  calculations. 
3.3.8  Functional  Block  Diagram  of  a  Sync- Recovery  Circuit 

The  block  diagram  of  a  sync  recovery  circuit  for  the  specific 
case  of  Example  1  in  Section  3.3  7  is  shown  in  Figure  3.37,  In  this 
case 

n  =  11,  k  =  13,  and  L  =  2047 
PA(x)  =  X11  +  X2  +  1 
PB(x)  =  x11  +  x6  +  X5  +  X  +  1 
and  Bp  =  B315  =  Bj0  +  B7  +  Bz  +  Bj  +  BQ. 

The  feedback  connections  for  PN  generator  A  are  omitted  in  the 
diagram;  but  the  feedback  connections  and  the  phase-shift  network 
connections  for  PH  generator  B  and  the  lower  buffer  register  are  redrawn 
and  shown  in  Figure  3.38. 

In  the  normal  operating  mode,  all  the  connections  on  PN  generator  B 
are  disconnected  and  the  generator  acts  as  a  buffer  register  receiving  a 
B-digit  every  13  bit-rate  clock  times.  At  the  instant  when  sync  loss 
is  detected,  the  circuit  enters  into  the  sync  recovery  mode  and  the  clock 
is  changed  to  the  "recovery  rate"  which  is  much  faster  than  the  normal 
bit  rate  clock.  The  recovery  clock  rate  will  be  discussed  in  a  later 
section.  For  convenience  in  the  discussion  we  set  the  time  as  t=0.  The 
recovery  circuit  functions  as  follows. 
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1.  At  t  =  0 

(a)  Switch  Sg  opens,  the  content  of  the  last  stage  of  the 
generator  is  designated  as  BQ  and  that  of  the  rest  of  the  stages  as 

Bj,  respectively  and  that  of  the  first  stage  as  B^Q.  The  digit  AQ  =  B^ 

is  ready  at  the  last  stage  and  Aj  =  Bp  appears  at  the  output  of  the 
phase-shift  network. 

(b)  Switch  SQ  closes  to  let  Ag  enter  PH  generator  A  or  a 
buffer  register  to  be  loaded  into  generator  A  later.  Switch  SQ  opens 
then  and  remains  open  throughout  the  recovery  mode. 

(c)  Switch  S^  closes  to  let  A^  =  Bp  enter  the  first  stage 

of  the  lower  buffer  register.  Switch  S-j  remains  closed  throughout  the 
recovery  mode.  At  each  clock  time  one  digit,  1  =  °>  1>  2,. ..,10 

and  j  =  1,  2,..., 9,  enters  the  lower  register. 

2.  At  t  =  tp,  tp  =  lOj,  j  =  1,  2,. ..,9,  the  lower  buffer  register 
is  filled  up  with  the  following  contents: 

Bjp+10’  bjp+9’---  Bjp+r  Bjp  =  V 

Generator  B  is  cleared  and  switch  Sp  closes  allowing  the  contents  of 
the  lower  buffer  register  to  be  loaded  into  generator  B. 

3.  At  t  =  t^  where  tft  =  tp  =  lOj,  switch  S^  closes  to  let 

Aj  =  Bjp  shift  into  generator  A.  Switch  S^  opens  until  j  increments. 

4.  At  t  =  tp,  tp  =  90,  A-jg  appears  at  the  output  of  the  phase- 

shift  network;  and  switch  Sp  closes  to  let  A10  enter  generator  A.  The  initial 
state  to  generator  A  is  filled  up  and  generator  A  begins  to  shift  until 
sync  is  caught  up. 

3.3.9  Estimation  of  Sync-Recovery  Time 

When  a  declaration  of  sync  loss  is  made,  let  ip  be  the  number 
of  bit  times  elapsed  after  the  most  recent  B-digit,  namely  B^  -j ,  entered 
generator  B  and  let  i^  be  the  corresponding  number  when  sync  is  recovered, 
as  illustrated  in  the  recovery  timing  diagram  in  Figure  3.39. 
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It  will  take  generator  A 

k(n-l)  +  1  +  iR  =  131  +  i 

shifts  to  catch  up  after  the  initial  state  has  been  determined  and  it 
will  take 

(n-1) (n-2) 

shifts  for  generator  B  to  determine  the  initial  state.  Thus  the 
combined  number  of  shifts  required  to  recover  the  sync  is 
k(n-l)  +  1  +  (n-1) (n-2)  +  iR 

If  recovery  is  required  within  no  more  than  one  bit  period  of  the 
original  bit-rate  clock,  then 

iR » v  1 

where  is  the  number  of  bit  times  elapsed  after  the  most  recent 
B-digit  entered  generator  B. 

0  £  i^  £  k-1, 
and 

1  <  iR<  k. 

Therefore,  the  total  number  of  shifts  required  by  both  generator  A 
and  generator  B  to  recovery  the  sync  is 
k(n-l)  +  1  +  (n-1) (n-2)  + 1  £  Number  of  shifts  <_  k(n-l)  +  1  +  (n-1) (n-2)  +  k. 
The  recovery  clock  should  be  approximately 

k(n-l)  +  1  +  (n-1) (n-2)  +  k  =  (k  +  n-3)  n+3 
times  faster  than  the  bit-rate  clock.  Thus,  if  R  is  the  bit  rate,  then 
the  recovery  clock  rate  is 

[n(k  +  n-3)  +  3]R. 

For  large  values  of  n  and  k,  the  recovery  clock  rate  is  approximately 
n(k  +  n-3)R. 
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For  the  example  we  were  discussing, 
n  =  11  and  k  =  13 

the  recovery  clock  is  234  times  faster  than  the  bit-rate  clock.  The 
bit  rate  of  the  Spectral  Data  Transfer  System  is  on  the  order  of 
kilobits  per  second  and  the  recovery  clock  rate  will  be  below  megabits 
per  second  which  is  well  within  the  state  of  art. 
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3.4  Preamble 

3.4.1  Purpose  of  the  Preamble 

In  the  Spectral  Data  Transfer  System,  the  receiver  must  establish 
a  coherent  carrier  reference  in  order  to  demodulate  the  received  signal. 

It  must  also  establish  a  bit  timing  reference  in  order  to  provide  sampling 
pulses  to  the  data  demodulator  and  to  clock  any  baseband  digital  processing. 
Furthermore,  frame  synchroni zation  must  also  be  established  in  order 
to  identify  correctly  the  self-sychronization  bits  for  the  scrambler  and 
the  scrambled  data. 

In  order  to  permit  the  receiver  to  perform  these  essential 
functions  efficiently  and  effectively,  each  transmission  from  the  transmitter 
begins  with  a  special  preamble.  This  preamble  is  designed  to  facilitate 
the  receiver's  synchronization  tasks:  carrier  synchroni zation,  bit 
synh-roni zation,  and  initial  frame  and  PU  Generator  synchronization. 

3. 4. 1.1  Carrier  Synchroni zation 

The  receiver  employs  a  loop  structure  to  track  the  noisy 
received  carrier  and  generate  a  coherent  local  reference  for  the  demodulation 
of  the  data.  The  exact  details  of  the  carrier-tracking  structure  will 
depend  upon  the  type  of  modulation  impressed  upon  the  carrier  by  the 
transmitter.  If  the  modulation  suppresses  the  carrier,  i.e.  there  is  no 
spectral  line  at  the  carrier  frequency,  then  a  structure  capable  of 
generating  a  line  related  to  the  carrier  frequency  must  be  employed, 
rather  than  a  simple  tracking  loop. 

The  two  structures  most  commonly  employed  for  carrier  tracking 
or  suppressed  carrier  tracking  are,  respectively,  the  phase-locked  loop 
{ PLL )  and  the  Costas  loop.  These  structures  and  their  tracking  performances 
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are  treated  thoroughly  in  the  literature  [8]-[18].  We  will  quote 
the  results  necessary  for  the  design  of  the  carrier  recovery  portion  of 
the  preamble. 

The  frequency  of  the  received  carrier  may  differ  from  the 
quiescent  frequency  of  the  VCO  in  the  receiver's  carrier  tracker  due  to 
effects  such  as  oscillator  drift  and  Doppler  shift.  Therefore,  the 
tracker  must  be  able  to  lock  in  the  presence  of  the  maximum  expected 
frequency  offset.  For  a  second  order  loop,  which  is  often  employed 
for  carrier  tracking,  with  a  loop  filter  specified  by 


F<s>  ■ K  t4t?  i3-1’ 

the  pull-in  range  Au  is  given  approximately  by  [8,  p.  364]: 

Aw  -  2%^TT-  ( 3-2) 

where  ion  is  the  natural  frequency  of  the  loop. 

The  acquisition  time  T  is  approximately  [8,  p.  364]: 

dCC} 

Tacq  “  - 3 - TT~ - 2 -  “  if  <VK*°  (3_3) 

q  4Cd)3  -  [(AU)\,Z/K]  -  (2«J/K)  2tu,3  n 

where  c  is  the  loop  damping  factor  and  K  is  the  amplitude  scale  factor 
of  F(s)  in  (3-1).  The  time  given  by  (3-3)  is  for  pull-in  to  a  pha*e 
error  of  tt/2  radians.  The  additional  time  Tg  to  settle  to  a  phase  error 
in  the  absence  of  a  frequency  offset  is  approximately  [8,  p.  365]: 


2B 


1  tn-2- 


T6 


(3-4) 


where  Bn  is  the  one-sided  noise  bandwidth  of  the  loop. 

The  preceding  discussion  has  not  considered  the  effects  of  noise 
on  the  acquisition  process.  The  lack  of  analytical  results  on  the 
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acquisition  time  of  phase-locked  loops  is  due  largely  to  the  inability 
to  solve  the  nonlinear  differential  equations  describing  the  loop. 
However,  some  results  obtained  by  computer  simulation  have  appeared  in 
the  literature  [15],  [17].  From  these  results  the  acquisition  time 
(for  90 %  probability  of  acquisition)  in  the  presence  of  noise-  with  no 
frqquency  offset,  can  be  approximated,  for  loop  SNR  ^  20  dB,  by 


„  10  2.99 

acq  f  B 

M  n  n 


(3-5) 


where  f  =  w  /2ir.  If  the  loop  SNR  decreases  to  10  dB,  then  T  „ 
n  n  acq 

approximately  doubles. 

The  design  of  the  preamble  proceeds  in  the  following  manner. 

In  order  to  maximize  the  carrier  power  (and  hence  loop  SNR)  available 

during  the  acquisition  process,  the  first  segment  of  the  preamble  will 

be  an  unmodulated  carrier.  The  receiver  designer  will  employ  (3-2)  as 

part  of  his  design  of  the  loop,  taking  into  account  the  maximum 

anticipated  frequency  offset  due  to  oscillator  instabilities  in  both  the 

transmitter  and  the  receiver,  Doppler  shifts,  etc.  Commonly  the  factor 

C  appearing  in  (3-2)  is  set  to  t  =  1//2  as  a  good  compromise  between 

speed  and  stability  of  the  loop  [9,  p.  54].  Once  (3-2)  and  other  design 

criteria  are  considered,  a  suitable  u  will  be  chosen  by  the  receiver 

n 

designer.  Then  (3-5)  and  the  results  reported  in  [15]  and  [17]  can  be 
used  to  determine  the  duration  of  the  unmodulated  carrier  which  forms  the 


first  segment  of  the  preamble. 

3. 4. 1.2  Bit  Synchronization 

Once  the  receiver  has  acquired  carrier  synchronization,  it  is 
able  to  recover  a  baseband  signal  from  the  modulated  incoming  signal.  The 
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next  step  is  to  synchronize  the  local  bit-rate  clock  in  frequency  and 
time  (epoch)  with  the  bit  transitions  of  the  received  signal.  This 
process  is  commonly  accomplished  by  a  tracking  loop  structure,  and,  thus, 
is  somewhat  similar  to  the  carrier  synchronization  process.  However, 
the  bit  synchronizer  locks  to  a  line,  related  to  the  bit  rate,  in  the 
baseband  spectrum,  whereas  the  carrier  tracker  operates  on  a  line  in 
the  RF  spectrum.  (Either  case  may  require  a  nonlinear  operation  to 
generate  a  suitable  line  to  track.)  Again,  the  subject  of  bit  syn¬ 
chronization  has  been  covered  extensively  in  the  literature  [8],  [10], 

['ll],  [13],  [15],  [19]  and  only  the  pertinent  results  will  be  mentioned 
here. 

Two  bit  synchronizer  structures  which  may  be  employed  are  the 
in-phase/mid-phase  tracker  shown  in  Figure  3.40  and  the  absolute-value 
early/late-gate  tracker  shown  in  Figure  3.41.  The  in-phase/mid-phase 
synchronizer  has  better  noise  tracking  performance  if  the  window 
(integration  interval)  is  small  (less  than  1/4  of  a  bit  period)  and  has 
a  shorter  pull-in  time,  but  the  early/late-gate  synchronizer  is  less 
sensitive  to  D.C.  offsets  and  is  simpler  to  implement  [8,  p.  445]. 

Transition-tracking  loops  used  for  bit  synchronization  can  be 
analyzed  exactly  as  with  phase-locked  loops  [13,  p.  458].  Thus  the  results 
cited  above  for  the  carrier  tracking  loop  may  be  applied  to  the  design 
of  the  preamble  for  bit  synchronization.  There  are,  however,  two 
differences  to  be  considered.  First,  of  course,  the  analysis  must  use 
the  parameters  of  the  bit  synchronizer's  loop.  Second,  the  preamble 
segment  for  bit-synchronization  acquisition  will  consist  of  a  pattern 
yielding  a  maximum  bit-transition  density.  For  an  NRZ-L  baseband  format 
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MID-PHASE 

INTEGRATOR 


Figure  3.40  IN-PHASE/MID-PHASE  BIT  SYNCHRONIZER 
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Figure  3.41  ABSOLUTE-VALUE  EARLY/LATE-GATE 
BIT  SYNCHRONIZER 
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the  pattern  would  be  alternating  O's  and  l's,  whereas  for  a  differentially 
encoded  (or  NRZ-M)  format  the  pattern  would  be  all-1 's.  The  pattern 
length  will  depend  upon  the  loop  SNR  and  the  loop  bandwidth  of  the 
bit  synchronizer,  but  a  few  tens  of  bits  typically  would  be  expected 
to  suffice. 

3. 4. 1.3  Unique  Word 

The  receiver  must  have  some  means  of  determining  when  the 
acquisition-aiding  preamble  concludes  and  data  transmission  begins,  in 
order  to  initialize  the  frame  synchronization  and  the  PN  Generator 
synchronization.  Therefore,  the  final  segment  of  the  preamble  is  a  unique 
word  which  can  be  detected  by  the  receiver  to  mark  the  epoch  at  which 
the  preamble  ends  and  data  transmission  begins.  In  order  to  insure  that 
the  epoch  determined  by  recognition  of  this  word  is  the  correct  epoch, 
the  word  should  have  a  high  autocorrelation  peak  when  properly  aligned 
and  low  correlation  when  misaligned. 

One  class  of  words  having  such  a  property  is  the  class  known 
as  Barker  cooes.  These  codes  consist  of  a  finite-length  word  of  ±1 
symbols  and  have  a  unique  property  that  the  autocorrelation  function 
C^,  defined  by 

N-k 

ck =  Z)  xi  Vk  <3_6) 

i  =  l 

where  N  is  the  word  length,  k  is  the  shift,  and  x-  is  a  code  symbol, 
satisfies  |C^|  _<  1  for  k  f  0. 

The  definition  of  the  autocorrelation  function  given  in 
(3-6)  yields  an  ideal  shape  (high  peak  for  k=  0,  uniform  low  value  for 
k  f  0)  only  if  adjacent  code  symbols  are  assumed  to  be  zero.  If  the 
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adjacent  data  bits  are  not  zero,  but  rather  are  random  ± 1 ' s  (as  must 


be  the  case  in  a  binary  system),  then  correlation  "sidelobes"  can 
arise  from  a  partial ' match  with  the  random  data.  The  longer  the  unique 
word,  the  less  is  the  probability  of  matching  the  random  data,  and  thus 
the  performance  is  improved.  A  list  of  all  known  Barker  codes  is  given 
in  Table  3.2,  from  which  we  see  that  the  longest  Barker  code  is  13  bits.* 
Codes  not  possessing  the  ideal  autocorrelation  properties  of  Barker 
codes,  but  of  greater  length,  may  discriminate  better  against  random  data. 
For  example,  a  100-bit  sequence  with  £  2,  k  =  1,  2,  ....  99,  would 
provide  better  discrimination  against  random  data  than  the  longest-known 
13-bit  Barker  code  [20]. 

Neuman  and  Hofman  [20]  conducted  an  exhaustive  computer  search 
for  sequences  wi  th  good  correlation  properties  for  lengths  up  to  24  bits. 
The  Neuman-Hofman  codes  are  listed  in  Table  3.3.  The  maximum  sidelobe 
of  the  autocorrelations  of  those  codes  is  <  2  for  all  codes  except 
the  first  24-bit  code  for  which  =  3.  The  magnitudes  of  the  auto¬ 
correlations,  which  are  of  interest  if  a  sign  ambiguity  exists  in  the 
received  data,  are  listed  in  Table  3.4  in  the  same  order  as  the  codes  are 
1  is  ted  in  Table  3.3. 

Massey  [21]  has  investigated  the  probability  of  erroneous 
synchronization  for  optimum  and  correlation  detection  of  Barker  and 
Neuman-Hofman  codes  imbedded  in  binary  data  transmitted  over  the  additive 
white  Gaussian  noise  channel.  His  results,  obtained  by  computer 
simulation,  are  summarized  in  Table  3.5. 

In  choosing  a  unique  word  for  use  in  the  spectral  data  transfer 
system,  we  desire  the  shortest  preamble  which  will  give  adequate 


♦Proof  exists  that  no  Barker  codes  exists  for  length  14  <_  N  <  6083, 
and  it  is  conjectured  that  Barker  codes  longer  than  N  =  T3  eiTist  [20]. 
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BARKER  CODES 


LENGTH 

CODE  SEQUENCE* 

1 

+1 

2 

+1  +1 

+1  -1 

3 

+1  +1  -1 

4 

+1  +1  +1  -1 

+1  +1  -1  +1 

5 

+1  +1  +1  -1  +1 

7 

+1  +1  +1  -1  -1  +1  -1 

11 

+]  +i  +]  _7  _]  _7  +7  -7  -7  +7  -7 

13 

+1  +1  +1  +1  +1  -1  -1  +1  +1  -1  +1  -1  +1 

♦The  complement  code,  the  reflected  code  (reverse  order),  and  the 
complement  reflected  code  also  are  Barker  codes. 


i ABLE  3.3 

NEUMAN-HOFMAN  CODES 


LENGTH 

CODE* 

7** 

+1 

+  1 

+1 

-1 

-1 

+  1 

-1 

8 

+1 

+  1 

+  1 

+1 

-1 

-1 

+1 

-1 

8 

+1 

+1 

+  1 

-1 

-1 

-1 

+  1 

-1 

9 

+1 

+  1 

+1 

+1 

-1 

-1 

-1 

+  1 

-1 

9 

+  1 

+  1 

-1 

-1 

-1 

-1 

-1 

+  1 

-1 

10 

+1 

+  1 

+1 

+1 

+1 

-1 

-1 

+  1 

-1 

+1 

10 

+1 

+  1 

+  1 

+1 

-1 

-1 

+1 

-1 

+1 

-1 

11 

+1 

+  1 

+1 

+1 

-1 

+  1 

+1 

+  1 

-1 

-1 

+1 

11** 

+1 

+1 

+  1 

-1 

-1 

-1 

+  1 

-1 

-1 

+  1 

-1 

12 

+  1 

+1 

+1 

-1 

-1 

-1 

-1 

+  1 

-1 

-1 

+1 

-1 

12 

+  1 

+  1 

-1 

-1 

+1 

+  1 

+  1 

+  1 

+1 

-1 

+  1 

-1 

13 

+1 

+  1 

+1 

+1 

+1 

+  1 

-1 

-1 

+  1 

+  1 

-1 

+  1 

-1 

13 

+1 

+1 

+  1 

+1 

+  1 

-1 

+  1 

-1 

-1 

+  1 

+1 

-1 

-1 

14 

+1 

+1 

+  1 

+1 

-1 

-1 

+  1 

+  1 

-1 

-1 

+  1 

-1 

+1 

-1 

14 

+1 

+  1 

-1 

-1 

+1 

+  1 

-1 

-1 

-1 

-1 

-1 

+  1 

-1 

+1 

15 

+1 

+1 

-1 

-1 

-1 

-1 

-1 

+  1 

+  1 

-1 

-1 

+1 

-1 

+1 

-1 

15 

+  1 

+1 

+  1 

+  1 

-1 

-1 

+  1 

+  1 

-1 

+  1 

+  1 

-1 

+  1 

-1 

+  1 

16 

+  1 

+  1 

+  1 

+  1 

+  1 

-1 

-1 

+1 

+1 

-1 

-1 

+  1 

-1 

+1 

-1 

-1 

16 

+  1 

+  1 

+  1 

+  1 

-1 

-1 

-1 

+  1 

-1 

-1 

-1 

+1 

-1 

-1 

+1 

-1 

17 

+  1 

+  1 

+1 

+1 

-1 

+  1 

-1 

-1 

+  1 

+  1 

-1 

-1 

-1 

+1 

-1 

+  1 

-1 

17 

+1 

+1 

+  1 

+1 

-1 

-1 

-1 

-1 

+  1 

-1 

-1 

+  1 

-1 

-1 

-1 

+1 

-1 

18 

+1 

+1 

+  1 

+  1 

-1 

+  1 

-1 

+  1 

-1 

-1 

+  1 

-1 

-1 

+  1 

+  1 

-1 

-1 

-1 

18 

+1 

+  1 

-1 

-1 

+  1 

+1 

-1 

-1 

-1 

-1 

-1 

+  1 

-1 

+  1 

+  1 

-1 

+  1 

-1 

19 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

+  1 

+1 

+  1 

-1 

+  1 

+  1 

+  1 

-1 

+  1 

+1 

-1 

+1 

19 

+1 

+  1 

+  1 

-1 

-1 

-1 

+  1 

-1 

-1 

-1 

+  1 

-1 

-1 

+  1 

-1 

-1 

+  1 

-1 

+  1 

20 

+1 

+1 

+1 

+1 

+  1 

-1 

+  1 

+  1 

-1 

-1 

+  1 

-1 

+  1 

-1 

+  1 

+  1 

-1 

-1 

-1 

+1 

20 

+1 

+1 

+  1 

-1 

+  1 

+  1 

+  1 

-1 

-1 

-1 

-1 

-1 

+  1 

+  1 

-1 

+  1 

-1 

-1 

+1 

-1 

21 

+1 

+  1 

+  1 

+  1 

+  1 

+  1 

-1 

+  1 

-1 

-1 

-1 

+  1 

-1 

+1 

+1 

-1 

-1 

-1 

+1 

+  1 

-1 

21 

+  1 

+  1 

-1 

-1 

+  1 

-1 

-1 

+  1 

+  1 

+  1 

+  1 

-1 

+  1 

+  1 

+  1 

+  1 

-1 

+  1 

-1 

+  1 

-1 

22 

+1 

+1 

1 

-1 

+1 

+  1 

+  1 

-1 

-1 

-1 

-1 

-1 

+  1 

+  1 

-1 

-1 

+  1 

-1 

-1 

+1 

-1 

+1 

23 

+1 

+1 

+  1 

+1 

+  1 

+  1 

-1 

+  1 

-1 

+  1 

-1 

-1 

+1 

+  1 

-1 

-1 

+1 

-1 

+1 

+  1 

-1 

-1 

-1 

23 

+  1 

+  1 

+  1 

+  1 

+  1 

+1 

-1 

-1 

-1 

-1 

+  1 

+  1 

-1 

-1 

+  1 

+  1 

-1 

+1 

+1 

-1 

+1 

-1 

+  1 

24 

+  1 

+1 

+  1 

-1 

-1 

-1 

-1 

-1 

-1 

+  1 

+  1 

-1 

+  1 

+1 

+1 

+  1 

-1 

-1 

+  1 

+1 

-1 

+1 

-1 

+1 

24 

+  1 

+1 

+1 

+  1 

+1 

-1 

-1 

-1 

+  1 

+  1 

-1 

-1 

-1 

+  1 

-1 

+1 

-1 

+  1 

-1 

-1 

+1 

-1 

-1 

+  1 

♦The  complement  code,  the  reflected  code,  and  the  complement  reflected  code  are  also 
Neuman-Hofman  codes. 

♦♦This  is  also  a  Barker  code. 


TABLE  3.4 

AUTOCORRELATION  S1DELOBES  OF 
NEUMAN-HOFMAN  CODES 


1 

1 


TABLE  3.5 

FALSE-SYNCHRONIZATION  PROBABILITY  FOR 
BARKER  AND  NEUMAN-HOFMAN  CODES 


FALSE-SYNCHRONIZATION  PROBABILITY  FOR 


CONDITIONS 

13-BIT 

BARKER 

CODE 

13-BIT 

NEUMAN-HOFMAN 

CODE 

7-BIT 

BARKER 

CODE 

NO  SIGN  AMBIGUITY 

OPTIMUM  DETECTOR 

E/Nq  =  1/2 

0.31 

0.28 

0.40 

1 

0.09 

0.07 

0.21 

2 

0.00 

0.00 

0.09 

CORRELATION  DETECTOR 

E/Nq  =  1/2 

0.42 

0.32 

0.45 

1 

0.19 

0.18 

0.32 

2 

0.08 

0.07 

0.22 

BPSK  WITH  SIGN  AMBIGUITY 

OPTIMUM  DETECTOR 

CVJ 

II 

o 

z 

LU 

0.39 

0.39 

0.63 

1 

0.14 

0.14 

0.37 

2 

0.00 

0.00 

0.21 

CORRELATION  DETECTOR 

E/Nq  =  1/2 

0.47 

0.49 

0.63 

1 

0.27 

0.24 

0.46 

2 

0.12 

0.13 

0.40 

performance  so  that  prime  power  is  not  wasted  on  unnecessary  transmission, 
as  well  as  maximizing  the  amount  of  useful  data  which  can  be  relayed 
through  the  satellite.  Looking  at  the  sync  performance  in  Table  3.5, 
we  can  conclude  that  a  13-bit  Barker  sequence  would  be  a  reasonable 
choice  for  the  unique  word. 

The  receiver  can  detect  the  unique  word  by  using  correlation 
techniques,  either  active  or  passive.  An  example  of  the  receiver's  unique 
word  detector  using  a  passive  correlator  is  shown  in  Figure  3.42.  When 
the  received  unique  word  is  in  the  register,  the  comparator  recognizes 
the  coincidence  of  the  received  unique  word  with  the  local  replica  of  the 
13-bit  Barker  code  and  outputs  a  signal.  This  signal  can  be  used  to 
enable  the  clock  to  the  receiver's  PN  generators,  to  enable  the  PN  sync 
recovery  circuits,  and  to  provide  a  reference  epoch  to  the  demultiplexer. 

3.4.2  Structure  of  the  Preamble 

From  the  discussion  in  Section  3.4.1,  we  arrive  at  the  structure  of 
the  preamble,  which  is  illustrated  in  Figure  3.43.  We  see  that  the  pre¬ 
amble  consists  of  3  segments,  sent  time-sequentially .  The  first  segment  is 
an  unmodulated  carrier  for  carrier-sync  recovery.  This  segment  lasts 
for  Tq  seconds.  It  is  followed  by  a  bit-sync  pattern  (alternating  1 1 s  and 
3's  or  all-l's,  depending  upon  the  baseband  coding  format,  as  discussed 
in  Section  3. 4. 1.2)  which  lasts  for  T^  seconds.  The  bit-sync  pattern  is 
followed  immediately  by  the  unique  word  for  frame  synchronization.  The 
unique  word  has  a  duration  of  T^  seconds,  after  which  time  the  data 
frames  are  sent. 

3.4.3  Generation  of  the  Preamble 

A  block  diagram  of  a  circuit  to  generate  the  preamble  is  shown 
in  Figure  3.44.  Although  this  circuit  and  the  subsequent  operational 
description  are  in  terms  of  a  specialized  hardware  implementation,  it 
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Figure  3.42  PASSIVE  CORRELATOR  FOR  UNIQUE-WORD  DETECTION 
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Figure  3.44  PREAMBLE  GENERATION 
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must  be  borne  in  mind  that  a  sof tware/fi rrnware  implementation  using  a 
microprocessor  is  a  viable  alternative. 

The  implementation  example  shown  in  Figure  3.44  assumes  the  use 
of  a  binary  DPSK  modulator.  If  other  modulations  with  higher-order 
alphabets  (e.g.  QPSK)  are  used,  then  the  details  would  differ  in  order 

to  generate  the  proper  symbol  patterns. 

The  principle  of  the  preamble  generator  is  to  select  sequentially 

inputs  to  the  modulator  which  create  the  three  segments  of  the  preamble. 
After  the  completion  of  the  preamble,  the  FEC-coder  outputs  are  connected 
to  the  modulator  so  that  the  FEC-coded  scrambled  data  may  be  transmitted. 
The  stepping  of  the  selector  from  one  source  to  another  is  controlled 
by  a  timing  circuit  in  accordance  with  the  preamble  design.  In  addition, 
the  bit-rate  clock  to  the  scrambler  is  inhibited  during  the  preamble 
transmission  so  that  the  proper  initial  load  is  present  when  scrambled 
data  is  first  sent  to  the  modulator  via  an  FEC  coder. 

The  generation  of  the  preamble  is  initiated  by  an  external 
START  signal  which  is  applied  to  a  2-bit  counter  and  the  preamble -ti mine 
generator.  (In  addition,  this  START  signal  would  turn  on  the  RF 
transmitter;  but  this  is  external  to  the  generation  of  the  preamble.) 

The  assertion  of  the  START  signal  causes  the  2-bit  counter  to  be  reset 
to  the  state  00  and  initiates  the  first  timing  interval  of  the 
preamble -timing  generator  which  provides  clock  pulses  to  the  2-bit 
counter. 

The  preamble -timing  generator  is  a  device  which  outputs  three 
pulses  in  responses  to  a  START  input  pulse.  The  timing  of  the  output 
pulses  is  shown  in  Figure  3.45.  The  design  is  such  that  the  leading  edge  of 
eacn  pulse  terminates  a  segment  of  the  preamble. 
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Figure  3.45  PREAMBLE  TIMING  SIGNAL 
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The  outputs  of  the  2-bit  counter  are  fed  to  the  control  inputs 
of  a  4-line-to-l-line  selector  and  also  to  three  count-decoding  circuits. 

f  These  decoding  circuits  are  designed  to  decode  counts  equal  to  01,  10, 

and  11,  respectively,  and  provide  a  signal  when  the  appropriate  count 
is  recognized.  (In  practice  the  three  decoders  could  be  merged  into 

i  one  2-line-to-4-line  decoder.) 

The  selector  operates  in  accordance  with  the  truth  table  shown 
in  Table  3.6.  The  output  of  the  selector  agrees  with  input  0,  1,  2,  or 
3  as  the  control  lines  B  and  A  are  00,  01,  10,  or  11,  respectively. 

Thus  the  selector  can  be  considered  as  a  four-position  switch.  The 
selector  controls  the  signal  going  to  the  modulator. 

Upon  the  receipt  of  a  START  pulse,  the  outputs  of  the  2-bit 
counter  are  reset  to  00  and  the  first  timing  interval  of  the  preamble¬ 
timing  generator  is  initiated.  The  outputs  of  the  three  count-detection 
circuits  all  remain  0  since  none  of  them  are  configured  to  detect 
count  00.  The  selector  selects  input  0,  which  is  hard-wired  to  a  logic- 
0  level.  This  produces  a  constant  logic-0  input  to  the  modulator,  and 

*  hence  yields  a  constant-phase  carrier,  which  constitutes  the  first 
segment  of  the  preamble. 

After  a  time  period  of  Tq  seconds,  the  preamble-timing  generator 

*  outputs  a  clock  pulse  to  the  2-bit  counter,  causing  the  counter  to 
count  up  to  state  01.  The  01  count  is  decoded  by  the  C0UNT=1  detector, 
which  outputs  a  signal  to  start  the  bit-sync  pattern  generator.  The 

'  presence  of  the  01  count  also  causes  the  selector  to  select  its  number  1 

input  which  is  connected  to  the  output  of  the  bit-sync  pattern  generator; 
hence  the  bit-sync  pattern  is  applied  to  the  modulator.  If  the 

*  modulator  employs  DPSK,  the  bit-sync  pattern  generator  is  designed  to  put 
out  all-1 's  data;  if  BP5K  is  employed,  then  the  bit-sync  pattern  generator 
is  designed  to  put  out  alternating  l's  and  0's. 
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TABLE  3.6 


TRUTH  TABLE  FOR  SELECTOR  GATE 


CONTROL 

B  A 

INPUTS 

0  12  3 

OUTPUT 

0 

0 

0 

X 

X 

X 

0 

0 

0 

1 

X 

X 

X 

1 

0 

1 

X 

0 

X 

X 

0 

0 

1 

X 

1 

X 

X 

1 

1 

0 

X 

X 

0 

X 

0 

1 

0 

X 

X 

1 

X 

1 

1 

1 

X 

X 

X 

0 

0 

1 

1 

X 

X 

X 

1 

1 

X  =  Don't  Care 
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After  a  time  period  o  T-j  seconds,  the  preamble-timing  generator 
outputs  a  clock  pulse  to  the  2-bit  counter,  causing  the  counter  to  count 
up  to  state  10.  The  10  count  is  decoded  by  the  C0UNT=2  detector,  which 
outputs  a  signal  to  start  the  one-shot  unique-word  generator.  This 
signal  also  serves  to  stop  the  bit-sync  pattern  generator.  Upon  being 
triggered  by  the  output  of  the  count-detection  logic,  the  one-shot 
unique-word  generator  outputs  the  chosen  unique  word  exactly  once, 
then  returns  to  a  quiescent  state.  The  presence  of  the  10  count  also 
causes  the  selector  to  select  its  number  2  input,  which  is  connected 
to  the  output  of  the  one-shot  unique-word  generator,  and  hence  sends 
the  unique  word  to  the  modulator. 

After  a  time  period  of  T^  seconds,  the  preamble-timing  generator 
outputs  a  clock  pulse  to  the  2-bit  counter,  causing  the  counter  to  count 
up  to  state  11.  The  11  count  is  decoded  by  the  C0UNT=3  detector,  which 
outputs  a  signal  to  enable  a  gate,  allowing  the  bit-rate  clock  to  be 
applied  to  the  data-scrambl ing  circuits.  This  same  signal  also  stops 
the  2-bit  counter  and  inhibits  further  counting  until  a  new  start  signal 
is  applied.  The  presence  of  the  11  count  also  causes  the  selector  to 
select  its  number  3  inputs,  which  is  connected  to  the  output  of  the 
FEC  coder  which  processes  the  scrambled  data,  thus  sending  the  FEC-coded 
scrambled  data  to  the  modulator. 

At  this  point  the  preamble  generation  has  been  completed  and  the 
system  placed  into  the  run  mode.  The  system  will  continue  to  operate 
until  shut  down  by  a  source  external  to  the  preamble  generator. 


» 
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4.0  TRANSMISSION  OF  DPSK  SIGNAL  THROUGH  SATELLITE  TRANSPONDER 
The  information  collected  by  the  spectral  detector  must  be 
transmitted  to  the  processor  for  analysis.  Typically,  a  satellite 
relay  communications  system  would  be  used  to  provide  remote  processing. 

A  modem  employing  differentially  encoded  phase-shift-keying  (DPSK) 
signal  transmission  with  differentially  coherent  demodulation  is  an 
appropriate  choice  when  the  circuit  simplicity  and  the  accompanying 
cost-effectiveness  is  the  over-riding  considerations  in  the  data 
communication  system  considered  here.  We  consider  such  a  DPSK  system 

operating  over  the  hard-limiting  satellite  channel. 

The  performance  of  the  communications  link  depends  upon  the  trans¬ 
mitter  and  antenna  on  the  platform,  the  satellite,  and  the  receiving 
shore-based  terminal.  In  order  to  give  a  quantitative  performance 
evaluation,  we  must  have  available  the  pertinent  specifications  of 
these  elements  of  the  link.  The  platform  containing  the  spectral  detector, 
could  be  a  buoy  in  a  real  situation. 

4.1  Buoy-to-Shore  Link 

4.1.1  Constraints  Imposed  by  the  Buoy  Environment 

A  small  buoy  does  not  permit  the  use  of  directive  antennas, 
since  it  is  an  unstable  platform  whose  location  is  time-varying  as  the 
buoy  drifts  in  the  ocean.  Therefore,  the  buoy's  antenna  is  considered 
to  be,  at  best,  a  hemispheric  pattern  with  at  most  3  dB  gain  over 
isotropic.  A  more  realistic  assumption  is  a  0-dB  antenna,  because  of  the 
varying  ground-plane  formed  by  the  ocean  surface,  feed  line  losses, 
antenna  imperfections,  etc. 

Since  the  buoy's  antenna  can  not  be  used  to  provide  gain,  the 
actual  power  out  of  the  final  amplifier,  less  feed  losses,  will  be  the 
EIRP.  Furthermore,  the  buoy  has  a  limited  prime  power  capability, 
which  places  limitations  on  the  feasible  transmitter  peak  power  of 
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around  50  Watts  maximum.  These  considerations  motivate  the  system  design 

to  use  the  lowest  practical  carrier  frequency,  in  order  to  minimize  the 
free-space  loss  [22,  p.  34-3]  L^$.  which  is  given  by 


Lfs  =  36.58  +  20  log1()f  +  20  log1()d  (4-1) 

where 

f  is  the  frequency  in  MHz,  and 
d  is  the  range  in  statutue  miles. 

The  choice  of  satellite  is  constrained  by  the  operating  frequency. 
Most  commercial  communications  satellites  operate  in  the  4  GHz/6  GHz  band, 
while  the  military  DSCS  satellites  operate  in  the  7  GHz/8  GHz  band.  At 
these  frequencies  the  free-space  losses  are  excessive  for  buoy  applications. 
This  leaves  only  FLTSAT  [23],  and  similar  satellites  such  as  GAPSAT 
[24,  pp.  57-63],  which  operates  at  the  UHF  (225-400  MHz)  band.  Typical 
parameters  for  a  UHF  satellite  are  given  in  Table  4.1.  Assuming  a  10° 
elevation  angle  to  the  satellite,  the  slant  range  would  be  25267  statute 
miles.  Using  (4-1),  we  then  find  the  uplink  free-space  loss  to  be  174.14 
dB  and  the  downlink  free-space  loss,  172.86  dB. 

The  ground  station  which  receives  the  signal  is  assumed  to  be 
of  low-to-moderate  cost;  A  typical  station  would  use  a  helix  antenna 
of.moderate  size  and  a  terminal  available  in  the  Navy's  inventory,  such 
as  an  AN/ARC-143B.  Such  a  terminal  and  antenna  combination  typically 
has  a  G/T  =  -14.8  dB/K. 

4.1.2  Link  Performance  Analysis 

The  performance  of  DPSK  transmission  over  a  hard-limiting  channel 
has  been  evaluated  for  the  practical  case  of  correlated  noise  and  SNR 
imbalance  at  the  phase  detector  [25].  The  results  of  [25]  express  the 

p  p 

bit  error  probability  as  a  function  of  uplink  SNR  (R^),  downlink  SNR  ( ) , 

p 

and  SNR  imbalance  at  the  phase  detector  (x  ).  Performance  curves  for 
2  2  2 

typical  values  of  R|j,  R^,  and  x  are  shown  in  Figures  4.1  and  4.2.  It 
was  found  [25]  that  the  noise  correlation  did  not  effect  the  results  when 
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TABLE  4.1 


Typical 
for  a 


Communications  System  Parameters 
UHF  Communications  Satellite 


Carrier  Frequency 
Uplink 
Downlink 

Satellite  Receiver  G/T 

Satellite  Transmitter  EIRP 

Orbit 

Transponder  Bandwidth 

Transponder  Type 

Access  Method 


299  MHz 
258  MHz 

-16.7  dB/K 

28  dBW 

Geosynchronous 

25  kHz 

Hard  limiting 

Single  carrier  per  transponder 
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the  a  priori  symbol  probabilities  are  equal.  Since  we  are  considering 
a  system  employing  a  data  scrambler,  this  assumption  is  valid  (the 
scrambler  tends  to  "randomize"  the  data). 

It  is  interesting  to  note  the  asymptotic  behavior  of  the  curves 

2 

in  Figures  4.1  and  4.2.  For  low  R^,  i.e.  the  downlink-limited  region 
2  2  2 

( Rd  <<:  Ru ) »  tfie  curves  for  a  constant  X  approach  the  same  limit  regardless 

2  2  2  2 
of  R^.  At  high  R^,  i.e.  the  uplink-limited  region  (R(j>>Ru)»  the  curves  for 

2  2 
a  constant  Ry  approach  a  constant  asymptote  regardless  of  the  imbalance  X  . 

2 

If  Rd  tends  to  infinity,  the  concept  of  an  SNR  imbalance  becomes  meaningless, 

and  the  performance  depends  only  on  uplink  SNR,  as  seen  in  Figure  4.2. 

2  2 

The  uplink  and  downlink  SNRs ,  Ru  and  R^,  in  Figures  4.1  and  4.2  are 
ratios  of  signal  power  to  noise  power.  Thus  the  receiver  bandwidths  of  the 
satellite  and  the  earth  station  implicitly  enter  the  calculations  through 
these  parameters.  For  the  uplink,  the  appropriate  bandwidth  is  the  noise 
bandwidth  of  the  transponder;  for  the  downlink,  the  receiver  bandwidth  is 
commonly  equal  to  the  bit  rate.  Under  these  assumptions,  we  find  that 
Ry  =  (E  -  6.22)  dB  where  E  is  the  EIRP  of  the  buoy's  transmitter  in  dBW  and 
Rd  =  (68.94  -  B)  where  B=10  log^Q(bit  rate).  If  the  power  output  of  the 

2 

buoy's  transmitter  is  in  the  range  of  1  Watt  to  50  Watts,  then  -6.22  f.  Ru  i  10.77 

2 

dB.  Also,  for  a  bit  rate  of  100  b/s,  Rd  =  48.94  dB  and  for  a  bit  rate  of 

2.4  kb/s,  Rd  =  35.14  dB.  It  is  quite  clear  that  the  uplink  will  be  the 

2 

controlling  factor.  Furthermore,  for  these  high  values  of  Rd  the  effects  of 
SNR  imbalance  at  the  phase  detector  will  be  negligible  and  the  curve  in 
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2 

Figure  4.2  labelled  =  »  applies;  thus  we  may  use  the  "classical" 
analysis  dicussed  below  to  a  good  approximation. 

When  the  effects  of  SNR  imbalance  on  the  link  performance  are 
negligible,  the  performance  of  a  single-carrier-per-transponder  link 
(i.e.  no  intermodulation  noise)  can  be  evaluated  by  use  of  the  equation 


[26] 


(4-2) 


where 


(no)jp  1 


s  the  uplink  carrier  power-to-noise  density  ratio  at  the 


output  of  the  hard  limiter.  If  (C/Nq)qqWN  is  large  compared  to  (C/Ng)^, 
then  (c/No^TOTAl=  (C/Ng)gp.  **  may  *5e  re^ated  to  the  carrier  power-to-noise 
density  ratio  at  the  input  to  the  limiter  [27]  by  the  factor  shown  in 
Figure  4.3.  Finally  the  total  C/Ng  and  the  bit  rate  are  used  to  calculate 
Efa/No  from  the  relation,  in  decibels. 


4  -10  >°9iolbit  rate)  (4'3> 

which  may  be  used  with  the  "classical"  DPSK  performance  curves  [13,  p.  303], 
[32,  p.  385]  to  assess  link  performance.  Using  (4-2),  (4-3)  and  Figure  4.3 
we  have  plotted  in  Figure  4.4  the  E^/Ng  at  the  shore  receiver  as  a  function 
of  the  buoy's  EIRP  for  data  rates  of  100  b/s,  1200  b/s,  and  2400  b/s. 

The  "classical"  DPSK  performance  in  terms  of  bit-error  probability 
is  given  by  [12,  p.  384] 

Pe  =  \  exp(-Eb/N0).  (4-4) 

By  applying  (4-4)  to  the  Eb/NQ  values  plotted  in  Figure  4.4,  we  can  plot 
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the  bit  error  probability  as  a  function  of  the  buoy's  EIRP,  as  is  done 
in  Figure  4.5. 

It  must  be  borne  in  mind  that  the  performance  given  by  (4-4)  and 
in  Figure  4.5  assumes  no  detuning  and  perfect  bit  synchronization.  The 
effects  of  Doppler  shift,  oscillator  drift,  noise  in  the  bit-rate  clock 
tracking  loop  (bit  synchronizer),  propagation  disturbances,  and  other 
deviations  of  a  practical  system  from  the  "ideal"  theoretical  system  must 
be  taken  into  account  in  addition  to  the  bit  error  probability  given  by 
(4-4)  or  Figure  4.5.  These  effects  require  that  a  "margin"  be  added  to 
the  required  buoy  transmitter  power  to  accommodate  fading  and  the  imperfections 
inherent  in  a  real  system.  This  margin  can  amount  to  several  decibels,  and 
will  be  a  function  of  the  required  link  performance. 

4.2  Shore-to-Buoy  Link 

We  now  turn  our  attention  to  the  link  from  the  shore  to  the  buoy. 

This  link  may  be  required  for  several  reasons,  such  as  a  need  to  command  the 
buoy  to  different  processing  modes. or  polling  in  a  TDMA  system.  We  assume 
this  link  will  be  supported  by  the  same  satellite  as  the  buoy-to-shore  link. 

The  shore  transmitter  is  assumed  to  be  similar  in  capability  to  an 
AN/ARC-143B  used  with  a  helix  antenna.  The  shore  station's  EIRP  is  assumed 
to  be  36  dBW.  The  satellite  parameters  are  as  given  in  Table  4.1.  The 
G/T  of  the  buoy  receiver  is  taken  as  a  parameter,  since  it  will  depend  upon 
the  antenna  structure  on  the  buoy  and  the  noise  figure  of  the  receiver’s 
input  stage. 

For  an  elevation  angle  of  10°,  the  uplink  and  downlink  free-space 
losses  were  found  previously  to  be  174.14  dB  and  172.86  dB,  respectively. 

The  uplink  C/Nq  is  73.76  dBHz.  For  a  25-kHz  transponder,  the  input  pr/NgB 
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to  the  transponder  is  29.78  dB,  and  Figure  4.3  shows  the  limiter  factor 
to  be  3  dB.  Thus  (C/Nq)[JP  in  (4-2)  is  76.76  dBHz. 

For  the  downlink  to  the  buoy,  (C/NQ)D0WN  =  (83.74  +  G/T)dBHz  where 
G/T  is  the  figure  of  merit  of  the  buoy's  receiving  system.  Since  a  simple 
omnidirectional  antenna  is  used  and  since  the  receiver  must  be  low  cost,  a 
reasonable  range  for  G/T  would  be  -40  dB/K  £  G/T  <_  -20  dB/K.  This  gives 
a  range  of  43.74  dBHz  <_  (c/Nq)Dqwn  i  63.74  dBHz  for  use  in  (4-2).  For  this  range 
of  (C/Nq)DqWN,  Rj  >  16.75  dB  for  bit  rates  of  500  bits/second  or  less. 

Again,  this  shows  that  the  "classical"  analysis  is  a  good  approximation 
to  the  link  performance.  Therefore  (4-2)  applies  to  the  shore- to-buoy 
link.  The  results  are  shown  in  Figure  4.6,  which  also  indicates  the 
link  margin  for  a  bit-error  probability  of  10"*\  Because  of  the  low  data 
rates  on  the  shore-to-buoy  link,  the  margin  is  large  and  the  link  will  not 
be  a  constraint  on  system  operations. 
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Figure  4.6  EXAMPLE  OF  TOTAL  COMMUNICATION-LINK  PERFORMANCE 
(SHORE-TO-BUOY) 
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5.0  FORWARD  ERROR  CONTROL  SCHEME 


As  discussed  in  Section  4,  the  command  link  (shore-to-buoy)  is  not 
a  constraint  on  the  system  performance  since  a  large  margin  in  signal-to- 
noise  energy  ratio  (Eb/Ng)  is  allowed  as  seen  in  Figure  4.6.  However,  as 
we  can  see  from  Figure  4.5,  the  data  link  (buoy-to-shore)  may  be  very  likely 
subject  to  the  power  or  bit-rate  constraint.  For  example,  in  order  to 
support  a  link  with  rate  100  b/s  at  a  bit-error  rate  of  10”^,  the  required 
buoy's  EIRP  is  only  0.3  watts,  while  for  a  bit  rate  of  2400  b/s,  it  is 
about  5  watts.  In  other  words,  reliable  data  transmission  is  highly 
dependent  upon  the  transmitter's  (buoy's)  power  and  the  bit  rate  used. 

Under  such  ci rcumstances ,  coding  techniques  can  play  an  important  role 
in  achieving  the  required  reliability  of  the  data  transmission  by  controlling 
channel  error  through  tradeoff  between  power  and  bandwidth. 

There  are  two  general  classes  of  error  control  techniques.  When 
the  channel  is  two  way  (not  necessarily  full  duplex),  the  transmitter  can 
store  the  data  until  a  verification  or  a  request  for  repeat  is  received 
for  a  data  block.  This  system  is  called  "Automatic  Request  for  Retransmission" 
(ARQ).  Itneedssome  redundancy  for  error  detection,  and  requires  storage 
at  the  transmitter.  The  satellite  channel  has  a  substantial  time  delay  on  the 
order  of  0.25  second.  This  makes  ARQ  costly  if  not  impractical. 

The  second  class  is  "Forward  Error  Correcting"  (FEC)  codes.  In  this 
system  a  redundancy  is  built  into  the  data  to  provide  the  decoder  at  the 
receiver  with  error-correcting  capability.  When  FEC  coding  is  used,  the 
signal-to-noise  ratio  (Eb/Ng)  required  to  achieve  a  given  performance  can 
be  reduced  through  error  correction,  by  as  much  as  5  or  6  dB  in  practice. 

Thus,  the  signal  power  required  to  sustain  a  given  error  rate  can  be  reduced 
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significantly.  This  reduction  in  required  received  power  level  is  termed 
"coding  gain." 

There  are  two  types  of  coding  techniques  most  commonly  used  in  FEC 
scheme:  block  coding  technique  and  convolutional  coding  technique.  In 
what  follows  we  will  briefly  describe  these  two  coding  techniques  and  show 
performance  (coding  gain)  curves  of  some  candidate  codes  which  are  applicable 
to  the  Spectral  Data  Transfer  System. 

5.1  Block  Codes 

A  block  encoder  stores  k  data  bits,  encodes  them  into  n  code-bits  (n>k), 
and  sends  them  to  the  modulator;  the  digits  in  a  block  of  n  code-bits  are 
independent  of  those  in  either  the  preceding  blocks  or  the  succeeding  blocks. 
Block  codes  have  rigid  structure  based  on  mathematical  theory,  particularly, 
the  theory  of  finite  field. 

The  concept  of  (Hamming)  distance  is  useful  in  discussing  the  error- 
correcting  capability  of  codes.  The  "(Hamming)  distance"  between  two 
codewords  is  defined  to  be  the  number  of  positions  in  which  the  words  differ. 
Thus,  minimum  (Hamming)  distance  is  a  measure  of  error-correcting  capability. 
For  example,  if  the  minimum  distance  of  a  code  is  d,  then  it  is  possible  to 
correct  all  error  patterns  of  t  or  fewer  errors  if  and  only  if 
d  >  2t  +  1 . 

The  code  rate  R  of  a  block  code  is  defined  to  be  the  ratio  of  the 
number  of  data  bits  k  to  the  block  length  n,  k/n,  and  thus  indicate  the 
efficiency  of  a  code. 

A  class  of  codes  known  as  the  "BCH  codes"  have  powerful  error- 
correcting  properties  and  known  decoding  algorithms,  which  are  conceptually 
complicated  but  relatively  simple  in  terms  of  decoding  time  and  required 
circuitry  in  implementation.  In  these  codes,  long  codes  can  be  constructed 


with  both  code  efficiency  (code  rate)  and  error-correcting  capability.  On 
the  other  hand,  short  block  codes  are  attractive  from  implementation  point 
of  view.  Consistent  with  the  requirement  of  simple  implementation  for 
the  Spectral  Data  Transfer  System,  we  shall  show  in  subsection  5.5  the 
performance  curves  of  two  relatively  short  codes,  (15,  7)  BCH  code  and 
(24,  12)  Golay  code,  together  with  those  of  convolutional  codes. 

5.2  Convolutional  Codes 

Convolutional  codes  differ  from  block  codes  in  that  the  digits 
generated  by  the  encoder  in  a  particular  time  unit  depend  not  only  on  the 
information  digits  within  that  time  unit,  but  also  on  the  information 
digits  within  a  previous  span  of  time  units.  At  any  time  unit,  a  block 
of  kQ  information  digits  is  fed  into  the  encoder,  and  a  block  of  n^  code 
digits  is  generated  at  the  output  of  the  encoder,  where  kg<n0.  The  n^-digit 
output  block  not  only  depends  on  the  kQ-digit  information  block  of  the  same 
time  unit  but  also  depends  on  the  previous  K-l  information  blocks.  The  ratio 
kg/nQ  is  the  code  rate  and  Kk^  is  the  constraint  lenqth. 

The  error-correcting  capability  of  a  convolutional  code  is  determined 
by  its  minimum  distance.  However,  unlike  block  codes,  convolutional  codes 
do  not  have  mathematical  structure;  there  is  no  analytic  method  to  construct 
convolutional  codes  with  maximum  minimum  distance.  Some  progress  has  been 
made  in  developing  algorithms  for  finding  convolutional  codes  of  moderate 
length  with  good  distance  properties  [28]  [29],  and  computer-generated  codes 
have  been  tabulated  by  several  researchers  [5,  p.  411]. 

There  are  three  well-known  decoding  techniques  for  convolutional 
codes:  threshold  decoding,  sequential  decoding,  and  Viterbi  decoding.  Thres¬ 
hold  decoding  is  an  extremely  simple  technique  applicable  to  many  short 
codes  correcting  a  few  errors,  and  easily  extendable  to  correct  bursts 
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of  errors.  Its  efficiency  diminishes  as  the  number  of  errors  to  be 

corrected  becomes  large.  Sequential  decoding  is  the  best  performing 

practical  technique  for  achieving  very  low  probability  of  error.  For 

a  number  of  reasons  such  as  buffer  size  requirements,  computation 

speed,  and  metric  sensitivity,  sequential  decoding  is  not  suited  for 

high-speed  operations  and  is  more  complex  than  the  other  decoding  scheme, 

notably  the  Viterbi  decoding.  Viterbi  decoding  is  particularly  desirable 

for  efficient  communication  where  low  probability  of  error  is  not 

required.  The  crossover  point  above  which  Viterbi  decoding  is  preferable 

to  sequential  decoding  occurs  at  values  of  probability  of  error  somewhere 
-3  -5 

between  10  and  10  [30]  depending  on  the  transmitted  data  rate.  As 

the  data  rate  increases,  the  probability  of  error  crossover  point  decreases. 

5.3  Location  of  the  FEC  Encoder 

Forward  error  correcting  coding  can  be  implemented  in  three  different 
positions  in  the  system  as  shown  in  Figure  5.1.  In  position  I,  the  data  is 
encoded  before  scrambled  while  in  position  II,  the  data  is  scrambled  first 
then  encoded;  in  both  cases,  the  Y-sequence  is  encoded  separately.  Hence 
in  both  position  I  and  position  II,  two  separate  Encoder/Decoder  pairs  are 
required.  In  position  III,  the  Y-sequence  and  the  scrambled  data  sequence 
are  combined  before  being  encoded;  a  single  Encoder/Decoder  is  required  but 
it  operates  at  double  speed.  The  proposed  system  operates  at  a  transmission 
rate  in  the  kbps  range,  well  within  the  start  of  the  art;  hence  from  equipment 
complexity  point  of  view,  position  III  is  a  better  choice. 

5.4  Selection  of  Error  Correcting  Codes 

We  have  reviewed  two  important  coding  techniques  which  are  applicable 
to  the  spectral  data  transfer  system.  The  purpose  of  using  error  correcting 
codes  in  the  data  transfer  system  is  to  achieve  an  acceptable  reliability  in 
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Figure  5.1  Location  of  the  FEC  Encoder 
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face  of  power  limitation,  not  to  enhance  the  ultimate  channel  capacity; 
most  short  codes  with  simple  decoding  procedure  could  be  proper  candidates. 
We  will  show  in  the  following  subsection  the  performance  curves  of  various 
short  codes  (block  codes  and  convolutional  codes)  for  which  encoders  and 
decoders  are  relatively  simple  to  implement. 

In  comparing  block  codes  with  convolutional  codes,  it  is  widely 
claimed  that  convolutional  codes  give  better  performance  on  the  space 
channel  than  block  codes  of  the  same  order  of  complexity  [30],  [31].  With 
regard  to  decoding  techniques,  Viterbi  decoding  is  a  proven  concept,  widely 
applied  in  satellite  communication  systems,  and  does  not  require  major 
technical  and  hardware  breakthroughs.  Sequential  decoding  may  be  ruled 
out  because  of  its  complexity. 

We  will  thus  consider  convolutional  codes  with  Viterbi  decoding 
as  the  candidate  scheme  for  the  Spectral  Data  Transfer  System,  and  will 
show  more  detailed  performance  of  this  scheme. 

5.5  Performance  Comparison  of  Various  Coding  Schemes 

We  now  present  the  performance  curves  of  several  short  codes  of 
interest.  Figure  5.2  shows  the  bit  error  rate  P(e)  vs.  E^/Nq  curves  for 
(15,  7)  BCH  code,  (24,  12)  extended  Golay  code,  and  rate-1/2  convolutional 
codes  of  constraint  length  Kkg=k=3  and  7  with  Viterbi  decoding.  The  curve 
for  the  case  of  no  coding  is  plotted  from  the  DPSK  performance  (4-4): 

P(e)  =  1  exp(-Eb/NQ). 

As  mentioned  earlier,  coding  gain  is  the  reduction  of  E^/Ng  acquired 
by  coding  in  achieving  a  given  P(e).  For  example,  at  P(e)  =  10”^,  we  obtain 
from  Figure  5.2  coding  gains  of  0.9  dB,  1.6  dB,  2.2. dB,  and  3.3  dB  for 


167 


BIT  ERROR  RATE  P(e) 


tiaiigtit 


IliiiilllllilllllSiliiaifSlIliii"" 


Irate  1/2  convolutional  code 

WITH  VITERBI  DECODING  (K=3)L 


RATE  1/2  CONVOLUTIONAL 
CODE  WITH  VITERBI  DECODING 
(  - 


BSsaSSSsisSsSssilP^^^ssif5Sa»sss&sSsBssBlsSlisI=^ 


|  |  -  fiiiipiniii 


;  :  !>  •  !  !  -t  ■ 

7.7  TT-- 

-6  1  :  I  !’  fflT 


4  5  6  7  8  9  10  11  12  13  14 

Eb/N0  (dB) 

!:jrp  BIT  ERROR  RATE  VS.  Eb/N0  FOR  VARIOUS  CODING  SCHEMES 
(HARD  DECISION  -  Q  =  2) 


168 


(15,  7)  BCH  code,  rate  1/2  convolutional  code  of  K=3,  (24,  12)  Golay  code. 


and  rate  1/2  convolutional  codes  of  K=7  respectively.* 

Convolutional  codes,  in  connection  with  the  simplicity  of  their 
encoders  and  decoders  for  short  constraint  lengths,  out  perform  block  codes 
of  the  same  order  of  complexity.  Note  that  the  performance  of  the  convolutional 
code  of  K=3  is  quite  close  to  that  of  (24,  12)  Golay  code  whose  decoding 
is  known  to  be  quite  messy.  For  this  reason,  we  consider  convolutional 
codes  of  short  constraint  lengths  with  Viterbi  decoding  as  the  candidate 
error  controlling  sheme  in  the  Spectral  Data  Transfer  System.  More 
discussion  on  the  performance  of  this  scheme  is  thus  in  order. 

The  coding  gain  of  a  convolutional  code  is  a  function  of  various 

parameters  such  as  number  of  quantization  levels  (Q),  constraint  length  (K),  code 

rate  (R),  and  decoding  delay  (or  path  history  length,  L).  Figures  5.3 

through  5.7  are  typical  examples  of  the  effects  of  these  parameters  on 

achievable  coding  gain  obtained  by  the  Viterbi  decoder.  Figures  5.3  and 

5-4  indicate  the  performance  improvements  that  can  be  realized  by  using  hard 

decision  (Q=2)  and  soft  decision  (several  quantization  levels)  as  compared  to 

the  uncoded  case  (no  coding)  for  K=3  and  5  respectively.  The  effect  of 

constraint  length  K  on  the  coding  gain  is  shown  in  Figures  5-4  and  5*5  for 

Q=2  and  8  respectively.  One  can  readily  see  that  there  is  an  increase  of 

-4 

about  0.3  dB  for  each  increment  in  K  at  a  bit  error  rate  of  10  .  In  order 

*The~i codTm,- gains  presented  here  are  obtained  from  the  classical  results 
available  in  the  literature  [32],  [33],  [34],  [35],  which  are  all  based  on 
binary  antipodal  signalling  or  PSK  modulation  with  coherent  demodulation. 

These  results  may  not  be  the  same  as  those  of  a  DPSK  system.  In  fact,  the 
performance  curves  of  PSK  and  DPSK  systems  for  the  uncoded  case  are  quite 
different.  However,  it  can  easily  be  shown  that  the  "coding  gains"  are 
approximately  the  same  for  both  cases  over  the  range  of  interest 
[P(e)  10-3  to  10“6].  Thus  we  may  use  the  results  of  the  PSK  system  directly 
on  the  DPSK  system.  The  net  effect  is  to  shift  all  the  curves  of  the  PSK 
case  by  a  certain  amount  of  Eb/Ng  due  to  the  difference  between 

Q(/2Eb/Ng)  where  Q(a)  =  — — exp(-x2/2)  dx  and  ^  exp(-Eb/Ng). 
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BIT  ERROR  RATE  VS.  E./Nn  FOR  RATE  1/2  VITERBI  DECODING 
HARD  QUANTIZED  RECEI vEDUDATA  (Q=2).  (PATH  HISTORY 
LENGTH  =32  BITS),  K=3,5,7 


BIT  ERROR  RATE  P(e) 


2  3  4  5  6  7  8  9  10  11  12  13 
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(dB) 
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Figure  5.6  BIT  ERROR  RATE  VS.  E^/Nq  FOR  RATE  1/2  VITERBI  DECODING. 

8-LEVEL  QUANTIZATION  (Q=8)  (PATH  HISTORY  LENGTH  =  32  BITS)  , 
K=4,  6,  8. 
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Figure  5-7  BIT  ERROR  RATE  VS.  Eb/NQ  FOR  RATE  1/3  VITERBI  DECODING, 

8-LEVEL  QUANTIZATION  (Q=8)  (PATH  HISTORY  LENGTH  =  32  BITS) 
K=4,  6,  8 


to  see  the  effect  of  code  rate  R,  we  present  in  Figure  5.7  the  performance 
curves  of  a  rate  1/3  code.  Comparing  Figure  5*6  with  Figure  5»7  we  can 
see  that  the  latter  offers  a  0.3  to  0.5  dB  improvement  over  the  former 
for  fixed  K  in  the  range  shown.  The  "improvements"  mentioned  here  are 
obtained,  of  course,  at  the  expense  of  some  other  quantity.  Note  that  rate- 
1/3  coding  requires  3  times  he  bandwidth  of  that  of  an  uncoded  system  while 
rate  1/2  requires  only  twice  the  original  uncoded  system  bandwidth.  The 
larger  constraint  length  means  higher  cost,  and  thus  the  systems  evaluator 
must  consider  several  factors,  knowing  the  individual  parameter  influence, 
when  making  the  recommendations .  It  seems  that  constraint  length  K=3  to  5 
may  suffice  for  the  Spectral  Data  Transfer  System. 

In  Figures  5.3  through  5.7,  the  path  history  length  (L)  is  32  bits. 

It  should  be  the  engineer’s  knowledge  that  a  value  of  L  of  4  or  5  times  the 
code  constraint  length  is  sufficient  for  negligible  degradation  from  optimum 
decoder  performance. 
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Appendix  A 

We  have  described  the  correlator  output  as 

y(T)  =  u(T)  +  n(T),  (A-l ) 

in  which  T  is  understood  to  be  the  interval  over  which  the  output 

is  observed.  Because  of  the  integrating  effect  of  the  lowpass  filter-- 

indeed,  one  model  of  the  lowpass  filter  which  has  been  used  is  a  T-second 

integrator--it  is  possible  to  assert  that  the  noise  term  n(T)  is  nearly 

Gaussian.  Since  n(T)  is  defined  to  have  a  zero  mean  value,  we  need  only 

its  variance  to  describe  its  statistical  behavior. 

It  matters  how  large  Oy(Hon^Tj)  is  compared  to  the  quantization 

level  q  .  The  natural  interoretation  of  q  assumes,  for  example,  that 
p  p 

o„  «  q  . 

y  u 

For  moderate  values  of  W^T  and  CNR,  it  can  be  shown  that 


2 

.  _qy _ 

(2ttWn)4 


2ttWnT 


where 


Y  = 


(2*Wn)‘ 


6  = 


2irW. 


max 


e  - 


W, 


m  s  L 


( A- 3) 


17fi 


N...  ..V  v  \ 


For  y  =  0 


62 

2/*r 

l  4fi2 

2  1 

2,WNT 

(CNR)2 

(2irlfNT)2 

CNR  (CNR)2 

(A-4) 


These  expressions  are  for  the  no-offset  configuration.  The  offset 
case  is  handled  by  replacing  6  by  j(l+lS)- 

The  required  WT  and  CNR  to  insure  that  o^<  ^  can  be  found  from 
these  expressions,  and  are  given  in  the  following  tables. 

Table  A-1  describes  the  quantization  of  the  FM  correlator  output  when 
no  offset  is  used.  For  an  input  resolution  whose  normalized  value  is 
1/20,  as  previously  mentioned  there  need  to  be  output  (uniform) 
quantization  levels  whose  normalized  values  are  1/800.  However,  Table 
A-1  shows  that  an  effective  nonuniform  quantization  scheme  can  be 
implemented  by  grouping  numbers  of  output  cells  to  correspond  to  single 
input  cells.  For  example,  20  output  cells  (41-60)  correspond  to  input 
cell  5. 

Because  the  effective  size  of  these  groups  of  output  cells  increases 

in  the  same  proportion  as  the  noise  o^,  the  mininum  WT  to  insure  that 

o  >6o  remains  constant, 
vi  y 

Fewer  cells  would  be  required  if  the  input  range  is  restricted.  If 
we  desired  to  limit  the  number  of  output  cells  to  512,  for  example,  we 
would  simply  modify  the  table  as  follows: 

Input  Cells  Input  Range  Nominal  Effective  Q  Output  Cells 

16  . 775< | 6 | < . 7996  0.7937  15/200  (30)  481-511 

17  (6 | >.7996  -  -  -  "OVERFLOW" 
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TABLE  A-1  QUANTIZATION  SCHEME  FOR  FM  CORRELATOR 
OUTPUT,  NO  OFFSET 


Input 

Cells 

Input 

Range 

Nominal 

l«l 

Effective  Qy, 
Output  cells  size 

Output 

Cells 

At  CNR=10, 
Minimum  WT 
for  qp>6oy 

0 

0< { 6 | < . 025 

0.00 

1/800  (1) 

0 

- 

1 

.025<  j  6 1 <.075 

0.05 

1/200  (4) 

1-4 

31 

2 

. 075< | 6 | < . 125 

0.10 

1/100  (8) 

5-12 

31 

3 

.125< (6 |< . 175 

0.15 

3/200  (12) 

13-24 

31 

4 

. 1 75< | 6 | < . 225 

0.20 

1/50  (16) 

25-40 

31 

5 

. 225  < 1 6 1 c- 275 

0.25 

1/40  (20) 

41-60 

31 

6 

.275  <  j  6 | <- 324 

0.30 

3/100  (24) 

61-84 

31 

7 

. 325< | 6 | <- 375 

0.35 

7/200  (28) 

85-112 

31 

8 

. 375<161<<425 

0.40 

1/25  (32) 

113-144 

31 

9 

. 425<|6l<-475 

0.45 

9/200  (36) 

145-180 

31 

10 

.475<|6l<-525 

0.50 

1/20  (40) 

181-220 

31 

11 

.525<! 6 1 < . 575 

0.55 

11/200  (44) 

221-264 

31 

12 

. 575< | 6 | <- 625 

0.60 

3/50  (48) 

265-312 

31 

13 

.625< ] 6 | <.675 

0.65 

13/200  (52) 

313-364 

31 

14 

. 675<161<<725 

0.70 

7/100  (56) 

365-420 

31 

15 

. 725< 1 6 1<<  775 

0.75 

3/40  (60) 

421-480 

31 

16 

.775<|6l<<825 

0.80 

2/25  (64) 

481-544 

31 

17 

<825<  j  6  j  <<875 

0.85 

17/200  (68) 

545-612 

31 

18 

<875<161<<925 

0.90 

9/100  (72) 

612-684 

31 

19 

. 925< | 6 | < . 975 

0.95 

19/200  (76) 

685-760 

31 

20 

<975<161<1<025 

1.00 

1/10  (80) 

761-840 

31 
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This  nice  compatibility  between  uniform  input  and  uniform  output 

quantizations  does  not  hold  when  an  offset  of  c  =  D _  is  used.  Table  A-2 

max 

details  a  compromise  that  results  from  designing  a  quantization  scheme 
to  involve  128  output  cells  by  limiting  the  range  over  which  quantization 
is  performed.  The  numbers  in  the  table  are  based  on  the  following 
considerations.  If  the  range  over  which  A/D  conversion  is  performed 
is  limited  to  100a*,  this  corresponds  to  an  input  range  of  [ 6 ( <a .  Now 
to  require  that  the  number  of  quantization  levels  be  a  power  of  2  is  to 
require  that 

range (vi  =  2a  _  ,  (A-5) 

TI^  2  1  (A5> 


(2m-l)QD 
2  +  (2m-l)QD 


(A-6) 


Thus  when  =  1/20  and  seven  bits  of  quantization  have  been  selected, 
the  resulting  value  of  a  is  127/167  and  Table  A-2  gives  the  quantization 

intervals.  The  implementation  of  this  scheme  proceeds  as  in  Figure  2.2  of 
the  text  in  which  it  is  shown  that  the  correlator  output  u  must  be  conditioned 

by  the  operations 

u  =  G^p-pmin)‘  (A-7) 

These  operations  are  performed  so  that  the  appropriate  range  of  u  matches 

the  input  range  of  the  quantizer (expressed  as  zero  to  full  scale  =  F.S.). 

It  should  be  noted  also  that  the  equivalent  input  quantization  cells 
shown  in  Table  A-2  are  not  uniform  in  size.  Therefore,  the  W^T  product 
required  for  insuring  small  correlator  output  noise  will  vary  from  cell  to 
cell  as  shown  in  the  table.  The  design  choice  here  would  be  the  maximum 
required  value  of  W.T, 


TABLE  A- 2  QUANTIZATION  SCHEME  FOR  FM  CORRELATOR  OUTPUT 
USING  OFFSET  AND  LIMITED  QUANTIZER  RANGE 


Input 

Cells 

Input  Range 

Effective  Qp, 
Nominal 6  Output  cell  size 

Output 

Cells 

At  CNR=10 
Minimun  WT 
for  qp>6ay 

-16 

6<-. 78694 

- 

- 

"OVERFLOW(-)!’ 

- 

-15 

-.78694<6<-. 73666 

-.76048 

1/167 

0 

118 

-14 

-.73666<5<-. 69455 

-.71483 

II 

1 

166 

-13 

-,69455<5<-. 62423 

-.65758 

2/167 

2,3 

62 

-12 

-.62423<o<-. 56514 

-.59361 

II 

4,5 

86 

-11 

-,56514«5<-. 51316 

-.53842 

II 

6,7 

110 

-10 

-.51316<6<- .46623 

-.48915 

II 

8,9 

134 

-9 

-.46623<6<-. 42309 

-.44424 

II 

10,11 

158 

-8 

-.42309<6<-. 38297 

-.40270 

ll 

12,13 

182 

-7 

-,38297<6<-. 32726 

-.35451 

3/167 

14-16 

96 

-6 

-.32726<6<-. 27582 

-.30107 

II 

17-19 

112 

-5 

-,27582<6<-. 22780 

-.25143 

II 

18-20 

128 

-4 

-.22780<6<-. 18260 

-.20488 

II 

21-23 

144 

-3 

-.18260<6<-. 12595 

-.15380 

4/167 

24-27 

92 

-2 

-.12595<6<-. 07276 

-.09897 

II 

28-31 

104 

-1 

-.07276<6<-. 02247 

-.04728 

II 

32-35 

116 

0 

.02247< 6  <.02537 

.00174 

II 

36-39 

128 

1 

.02537< 6  <.07107 

.04847 

II 

40-43 

140 

2 

.07107<  6  <.12559 

.09867 

II 

44-48 

100 

3 

.  12559< 6  <.17759 

.15188 

II 

49-53 

109 

4 

.  17759< 6  <.22738 

.20274 

II 

54-58 

119 

5 

.22738<  6  <.27524 

.25154 

II 

59-63 

128 

6 

.27524< 6  <.32136 

.29850 

II 

64-68 

138 

7 

. 32 136< 6  <.37467 

. .34828 

6/167 

69-74 

104 

8 

. 37467< 6  <.42598 

.40056 

II 

75-80 

112 

9 

.42598< 6  <.47551 

.45096 

II 

81-86 

120 

10 

. 4 755 1< 6  <.52343 

.49966 

II 

87-92 

128 

11 

.52343<  <5  <.56989 

.54684 

II 

93-98 

136 

12 

.56989<  6  <.62241 

.59637 

7/167 

99-105 

107 

13 

.62241< 6  <.67329 

.64805 

II 

106-112 

114 

14 

.67329< 6  <.72266 

.69815 

II 

113-119 

121 

15 

. 72266< 6  <.76388 

.74339 

6/167 

120-127 

172 

16 

6>. 76388 

' 

"OVERFLOW^)" 

~ 
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